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ABSTTE^ACT 

This  thesis  consists  of  a  three-dimensional  numerical 
analysis  of  the  Jet  Aircraft  Hush  House  located  at  Naval  Air 
Station  Jacksonville,  Florida.  Utilizing  the  PHOENICS  Code 
allows  for  the  determination  of  the  aerothermal 
characteristics  including  velocity,  pressure,  enthalpy, 
turbulent  kinetic  energy  and  the  dissipation  rate  of  turbulent 
kinetic  energy  in  the  facility  during  testing  of  the  U.S. 
Navy's  F-4  (Phantom  II)  J-79-GE-8  gas  turbine  engine  with 
afterburner.  How  and  by  what  method  PHOENICS  arrives  at  this 
solution  is  discussed.  Of  greatest  importance  is  the  resulting 
"behavior"  of  the  aerothermal  system. 

Problems  encountered  using  the  PHOENICS  Code,  resulting 
numerical  solutions  to  the  particular  facility,  comparison  to 
actual  test  data  and  recommendations  for  further  applications 
of  the  PHOENICS  Code  are  presented. 
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I  -  XNTRODUCTION 

New  and  overhauled  aircraft  gas  turbine  engines  undergo 
stringent  performance  testing  before  returning  to  operational 
use.  To  allow  for  safe  and  effective  monitoring  of  the  engine 
and  aircraft,  gas  turbine  Test  Cell  and  Hush  House  facilities 
are  utilized  to  perform  monitoring  of  key  temperatures, 
turbine  speeds  and  flow  rates  associated  with  the  engine.  The 
test  facility  additionally  provides  distortion-free  air  flow 
to  the  engine  inlet  (eliminating  rain  and  wind  effects),  for 
lower  exhaust  temperatures  exiting  from  the  facility  and 
reduced  noise  pollution  to  the  surrounding  area.  Furthermore, 
these  facilities  allow  for  final  engine  "trim"  settings  to 
ensure  that  maximum  power  is  available. 

Currently  the  test  facilities  designed  and  built  for  the 
Navy  vary  widely.  The  various  facilities  include: 

1.  The  Run-up  pad--where  a  tied  down  aircraft  conducts 
power  checks  in  an  open  air  environment. 

2.  Sound  suppressor--located  at  the  Run-up  pad  to  reduce 
noise . 

3.  Open  Test  Stand  (OTS ) --located  outdoors  that  supports  an 
aircraft  engine  and  has  portable  instrumentation 
attached. 

4.  Hush  House--a  totally  enclosed  facility  allowing  around 
the  clock  testing  of  aircraft. 

5.  Test  Cell  (TO --an  enclosed  facility  allowing  testing  of 

various  engine  sizes  out  of  the  airframe. 


Design  engineers  have  used  experience  and  common  sense  to 
build  these  facilities  with  varying  degrees  of  success. 
Aerothermal  monitoring  supplements  the  design  of  the 
structures . 

With  the  tremendous  technological  advances  made  in  recent 
years,  much  higher  temperatures  and  thrusts  are  available 
compared  to  similar  size  and  weight  engines  of  just  a  few 
years  ago.  Considering  that  nearly  67%  of  the  Navy  Test  Cells 
and  Hush  Houses  are  at  least  25  years  old  it  is  evident  that 
a  new  generation  of  standardized  facilities  meeting  federal, 
state  and  local  noise,  pollution,  security  and  safety 
standards  is  required.  [Ref.l:p.3] 

The  Navy  Environmental  Health  Center  cited  total  claims 
against  the  government  in  1981  at  $114  million  due  to  physical 
injury,  hearing,  or  respiratory  damage  to  military  and 
civilian  personnel.  Additionally,  the  1981  ruling  by  the 
federal  court  system  that  military  gas  turbine  test  facilities 
must  comply  with  Section  233  of  the  Clean  Air  Act  accelerated 
the  need  for  modern  test  facilities  meeting  the  above 
requirements  plus  those  projected  for  the  future.  [Ref.l:p.3] 

The  goal  of  the  Navy  is  to  provide  for  standardized 
aircraft  and  engine  test  facilities  that  meet  both  safety  and 
environmental  concerns  while  providing  the  latest  technology 
to   monitor   the   engines.   To   accomplish   this,   complete 


development  of  test  facilities  with  reliable  operation  and 
maintenance  throughout  the  facility  life  cycle  is  necessary. 
[Ref .l:p.l] 

When  contemplating  the  inherent  design  criteria  for  the 
new  facilities,  problems  such  as:  turbulence,  non-concentric 
flow  patterns ,  multiple  jet  engine  aircraft  (including  varying 
temperature  and  thrust  arrangements),  supersonic  velocities, 
extreme  temperature  and  pressure  gradients,  and  awkward 
geometries  must  be  appraised.  [Ref. 5] 

Enigmas  associated  with  the  design  of  the  structure  may 
not  be  obvious.  Complete  information  dealing  with  the 
aerothermal  properties  is  a  must  in  basing  decisions  on  new 
multi-million  dollar  expenditures.  Therefore,  to  minimize  the 
possibility  of  an  erroneous  design,  the  PHOENICS  Code  is 
employed  first  on  documented  Test  Cell  and  Hush  House 
structures  having  actual  field  test  data.  After  verifying  the 
accuracy  of  PHOENICS  with  the  test  data,  the  Code  eventually 
could  be  used  at  a  Navy  installation  (such  as  the  Naval  Civil 
Engineering  Laboratory)  to  design  future  test  facilities 
providing  accurate  modeling  data. 

In  this  instance,  the  performance  of  a  Navy  J-79-GE-8  gas 
turbine  engine  located  at  Naval  Air  Station  Jacksonville, 
Florida  Hush  House  facility  is  investigated.  The  data  is  used 
to  determine  the  behavior  of  the  existing  structure  and  will 


be   employed   to   assist   engineers   in   the   research   and 
development  (R&D)  phase  of  the  facility  design  project. 


II-     PRC^T  .KM    FORMLTLATICasr 

A-     HUSH    HOUSE    CHARACl'KRIZATIOISI 

The  Hush  House  under  investigation  is  an  "Air  Force"  type 
located  at  Naval  Air  Station  Jacksonville,  Florida.  In  1980 
the  Air  Force  purchased  25  and  later  54  demountable  Hush 
Houses  with  Test  Cell  capability.  This  standardization 
accompanied  with  such  a  large  scale  purchase  and  the 
flexibility  of  aircraft  or  engine  test  capabilities  allowed 
for  a  cost  effective  combined  design.  Because  of  the  large 
scale  purchase  the  cost  per  Hush  House  was  a  relatively  low 
$2.2  million.  The  only  major  modification  performed  by  the 
Navy  to  the  Air  Force  design  was  inclusion  of  an  additional 
AFFF  fire  fighting  system.  The  Hush  House  design  does  not 
address  air  pollution  standards.  Some  type  of  air  pollution 
abatement  system  will  undoubtably  be  required  in  newer 
designs.  [Ref.l:p63 

The  general  dimensions  of  the  Hush  House  is  28.7  meters 
(94  feet)  wide  by  25.5  meters  (83.7  feet)  long  with  the  oval 
shaped  "augmenter  tube"  extending  an  additional  31.3  meters 
(103  feet)  from  the  rear  of  the  Hush  House  (See  Figures  2.1 
and  2.2).  The  aircraft  is  placed  so  that  the  engine  exhaust 
is  approximately  4.6  meters  (15  feet)  from  the  augmenter  tube 
entrance.   The  most   common  construction  material   for   the 
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structure  is  steel.  The  design  is  geometrically  simple 
allowing  for  mass  production.  The  air  flow  in  the  last  two 
inlets  nearest  the  augmenter  tube  provides  both  engine  intake 
and  cooling  air  for  the  surrounding  engine  area.  The  final  7 
meters  (23  feet)  of  augmenter  tube  is  angled  at  approximately 
45  degrees  running  from  the  floor  upward  to  provide  for  a 
smoother  directional  transition  of  gases.  Acoustic  damping 
materials  line  the  augmenter  tube,  inlets  and  aircraft  bay 
walls  to  reduce  exhaust  noise. 

The  Hush  House  is  used  in  support  of  the  Naval  Air  Rework 
Facility  (NARF)  located  at  the  Naval  Air  Station.  Although  the 
Hush  House  is  presently  utilized  for  engines  fully  mounted  in 
an  airframe,  it  is  anticipated  that  engine  "Test  Cell"  testing 
(where  an  engine  is  mounted  on  a  test  stand)  will  be  approved 
in  the  future.  The  facility  is  normally  in  operation  from  0700 
to  2400  Monday  through  Friday  with  an  average  of  1  to  5 
aircraft  being  tested  per  day.  One  of  the  most  common  aircraft 
tested  is  the  F/A-18  Hornet  powered  by  two  General  Electric 
F-404  gas  turbine  engines.  Typically,  after  completion  of 
overhaul  of  the  aircraft  and  engines,  the  aircraft  will  be 
towed  to  the  facility  and  backed  into  the  Hush  House.  The  hugh 
aircraft  inlet  doors  to  the  facility  will  then  be  closed  to 
provide  an  enclosed  area  free  of  foreign  matter,  outside  wind, 
rain,  and  other  environmental  and  personnel  concerns. 


with  the  aircraft  fully  enclosed  and  securely  tied  down, 
a  normal  crew  of  three  mechanics  will  perform  final  engine 
"trim".  One  person  is  located  in  the  Hush  House  Control  Booth 
to  ensure  manning  of  the  fire  fighting  systems.  The  second 
will  physically  adjust  the  fuel  control  mechanism  on  the 
engine  via  radio  signals  received  from  the  cockpit  where  the 
third  crewman  is.  The  engines  are  controlled  during  testing 
by  the  third  crewman  in  the  aircraft  who  is  additionally 
responsible  for  monitoring  the  various  fuel,  oil,  %  power, 
hydraulics,  and  other  gauges.  Testing  normally  lasts  for  a 
period  of  two  hours.  Only  one  engine  is  tested  at  a  time  with 
the  other  normally  at  idle.  Air  intake  and  cooling  is  provided 
by  five  inlet  areas  located  on  either  side  of  the  Hush  House 
(See  Figures  3.2  and  3.3).  Wire  screening,  air  diffuser  and 
noise  abatement  equipment  is  present  at  the  inlets.  Installed 
fire  fighting  and  wash-down  equipment  is  present  in  the 
aircraft  bay  area  with  the  air  inlets  capable  of  being 
autom.atically  closed. 

B-   ISSUES  CONCEFUNTENG  THE  HUSH  HOUSE 

The  Hush  House  is  currently  performing  as  designed  with 
the  below  exceptions.  Investigation  indicates  some  cracking 
of  the  welds  which  fuse  together  the  stainless  steel 
"corrugated"  sheet  metal  overlaying  the  sound  absorption 

material  located  in  the  augmenter  tube.  Also,  actuation  of 


heat  sensors  in  the  overhead  of  the  aircraft  bay  sometimes 
occurs  when  the  engines  are  idled  for  a  period  of  time. 
Apparently  some  of  the  exhaust  gases  travel  upward  when  low 
mass  flow  rates  are  present  instead  of  going  down  the  31  plus 
meters  (103  feet)  of  augmenter  tube. 

In  general,  personnel  at  the  Naval  Air  Station  were 
pleased  with  the  overall  operation  of  the  Hush  House.  Although 
vibration  of  the  entire  physical  structure  has  been  noted  in 
this  facility  as  in  other  designs,  it  does  not  appear  to  cause 
any  structural  problers  such  as  the  requirement  for  inspection 
and  retightening  of  bolts  as  occurs  in  other  Navy  Hush  House 
designs . 

Physical  test  data  is  available  for  an  identical  Air  Force 
facility  for  the  F-4  Phantom  II  and  is  included  in  parts  of 
this  report  for  comparison  with  results  from  the  PHOENICS 
Code.  The  test  data  provides  an  important  comparison  for  this 
thesis  to  verify  the  computer  model.  [Ref.2] 

A  physical  inspection  of  an  identical  Hush  House  to  Naval 
Air  Station  Jacksonville's  was  under-taken  in  June  1988  at 
Dobbins  Air  Force  Base  in  Atlanta  Georgia.  Discussions  with 
the  various  Air  Force  personnel  corroborated  the  previously 
mentioned  problems  with  the  exception  of  the  heat  sensor 
activation  problem.  The  Dobbins  Air  Force  Base  Hush  House  was 


utilized   for   both  Test   Cell   and   Hush   House   operations 
involving  the  F-15  (EAGLE)  fighter. 


Exhaust 


Augmenter  Tube 


Air  Inlets 


Figure  2.1   Actual  Hush  House  (Top  View  Looking  Down) 


10 


Exhaust 


Augmenter  Tube 


Aircraft  Entrance  Doors 
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Figure  2.2   Actual  Hush  House  (Side  View! 
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HI.  THE  PHOEtvnCS  CODE 

A-  AN  EXPLANATION 

PHOENICS,   which  stands   for   "Parabolic  Hyperbolic  Or 

Elliptic  Numerical  Integration  Code  Series",  is  a  computer 
code  capable  of  simulating  fluid  flow  (including  changes  of 
state),  heat  transfer  and  chemical  reactions.  The  code  is 
copyrighted  by  CHAM  (Concentration  Heat  and  Momentum) 
Limited  which  was  founded  by  D.  Brian  Spalding  in  1972. 
[Ref . 3:p.l.l] 

As  shown  in  Figure  3.1,  the  PHOENICS  Code  is  structured 
into  three  major  components  consisting  of  "EARTH", 
"SATELLITE",  and  "GROUND  STATION".  The  "EARTH"  program  may  be 
considered  as  a  "black  box"  which  contains  the  computational 
procedures  for  all  applications.  The  "SATELLITE"  program 
consists  of  several  data-input  methods  which  define  the 
boundaries  and  scope  of  the  problem.  Finally,  the  "GROUND 
STATION"  program  has  various  subroutines  containing  built- 
in  physical  models  which  may  or  may  not  be  selected  depending 
on  the  program  application. 

By  default  "EARTH"  can  solve  for  up  to  25  dependent 
variables.  In  the  Hush  House  application  these  include 
pressure,  all  three  dimensional  velocity  components, 
enthalpy  (temperature),   turbulent  kinetic  energy,  and  the 
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turbulent  kinetic  energy  dissipation  rate.  The    dependent 
variables   can  be  enlarged  if  desired. 


PHOEMICS  input  library 


Figure  3.1   Satellite,  Earth  and  Ground 

The  PHOENICS  Code  is  exceptionally  flexible  allowing  for 
varying  sizes  and  scales  of  structures  from  "micro"  to 
"macro".  A  "final"  solution  is  highly  dependent  on  the 
accuracy  required  by  the  user  as  well  as  the  grid  sizing  used. 
The  higher  the  accuracy  and/or  number  of  "cells"  desired,  the 
greater  the  computational  effort  required  for  a  solution.  In 
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this  investigation,  a  "grid  independent  solution",  where  a 
sequence  of  "finer"  grids  would  be  used,  was  not  possible  due 
to  time  limitations. 

In  order  to  solve  a  given  problem  the  user  must  create  a 
"Ql"  file.  The  "Ql"  file  describes  all  aspects  of  the  system 
being  investigated.  This  file  is  partitioned  into  several 
groups  allowing  the  user  to  specify  various  aspects  of  the 
physical  system.  Examples  are  structure  size  and  shape. 
Additionally,  other  groups  are  capable  of  providing  the 
necessary  grid  spacing,  boundary  conditions,  engine  mass 
flow  rate,  engine  exhaust  temperature  and  printing  options. 
The  "Ql"  file  is  read  by  the  "SATELLITE"  program  which  is 
recognized  and  processed  by  "EARTH".  The  "EARTH"  program 
solves  discretized  conservation  equations.  A  typical 
conservation  equation  (before  discretization)  which  "EARTH" 
solves  can  be  expressed  in  the  form: 


5 

--    {rp<J>)    +    div{rpV<f)    -    rf*    grad   <l>)    =    rS» 
5t 


where : 

r  =  phase  volume  fraction 

<f>  £  dependent  variable 

r#  =  exchange  coefficient  (laminar  or  turbulent) 

Si  =  source  or  sink  term 
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p   =   density 

V   2   velocity  vector 

Here  $  is  the  dependent  variable  currently  in  question.  For 
example,  using  the  conservation  of  energy  equation,  *  would 
be  enthalpy  and  r#  is  the  diffusion  constant.  [Ref . 3 :p. 2 . 4] 

The  values  computed  and  displayed  for  each  cell  are  the 
solutions  of  the  algebraic  equations  relevant  to  the  laws  of 
physics  for  that  cell.  The  laws  of  physics  are  those  of 
conservation  of  mass,  energy  and  momentum  and  any  other 
dependent  variable.  The  general  form  of  the  algebraic 
linearized  equations  resulting  from  a  discretization  of  the 
previous  field  equation  is: 

aE<J>E    +    aw<{>K    +    aN<i>N    +    as4>s    +    au^a    +    aL<i>L    +    aT<i>T    +    S 
$,    =    

aE    +    aw    +    aN    +    as    +    an    +    at    +    ai    +    ap 

where  <^i.  is  the  value  of  the  dependent  variable  at  a 
particular  cell  and  subscripts  E,  W,  N,  S,  H  and  L  are 
neighboring  locations.  The  values  of  4>  at  these  locations  are 
known  from  a  previous  iteration.  Quantities  S  and  ap  express 
the  influence  of  a  source  of  the  entity  *.  The  a's  are 
influence  coefficients,  temporarily  treated  as  constants,  with 
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the  above  subscripts  expressing  interaction  between  neighbor 
cells.  [Ref .3:p.2.4] 

Solving  the  multitude  of  algebraic  equations  is 
complicated  and  time  consuming.  An  iterative  manner  is  used 
until  the  imbalance  between  the  left  and  right  sides  of  all 
the  above  equations  are  negligible.  Several  schemes  can  be 
used  to  perform  the  "sweeps"  (where  a  sweep  solves  for  the 
dependent  variables  for  each  cell  in  the  grid  one  time) 
necessary  for  the  iteration  process.  The  results  will  be 
identical  although  the  time  to  complete  the  task  will  vary 
with  the  iteration  scheme.  [Ref . 3 :p . 2 . 5] 

-  '  There  are  several  ways  to  determine  whether  the  problem 
has  "converged"  to  a  solution.  Results  are  monitored  to 
determine  whether  values  of  dependent  variables  are  remaining 
constant  or  fluctuating.  The  "error"  (residuals)  is  tabulated 
for  each  sweep  at  the  end  of  the  computer  printout .  By 
observing  "residuals"  from  each  sweep  a  pattern  of  decreasing 
("converging")  values  should  be  present.  If  the  "residuals" 
are  increasing  ("diverging")  adjustments  of  several  constants 
which  in-turn  apply  a  dampening  effect  upon  the  variables  in 
question  may  be  required.  The  printout  of  successive  residual 
values  demonstrate  convergence  by  a  reduction  in  the 
residuals.  That  is,  a  visually  downward  sloping  line  on  the 
accompanying    graph    of    residuals    shows    convergence. 
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Conservation  balances  on  each  dependent  variable  can  be 
verified  using  "net  source"  values  which  are  tabulated  at  the 
end  of  each  run.  The  mass  flow  rate  (mass  flow  rate  =  density 
*  area  *  velocity)  is  additionally  calculated,  in  this 
instance,  for  the  five  inlets  and  one  outlet  confirming  that 
no  hidden,  missing,  or  improper  boundaries  are  present.  The 
problem  has  "converged"  when  an  acceptable  level  of  error 
remains . 

Successful  "convergence"  of  the  PHOENICS  program  is 
dependent  upon  experience  and  an  understanding  of  the  results. 
If  the  residuals  of  a  dependent  variable  show  divergence,  the 
user  can  manipulate  the  rate  of  change  of  values  of  that 
variable  by  using  the  RELAX  command  to  return  the  variable 
back  to  the  convergence  path.  Relaxation  is  placed  on  a 
variable  to  control  the  variables  maximum  change  froni  each 
sweep.  This  minimizes  the  possibility  of  "over-correction" 
or  "blowing  up"  of  the  values.  Smaller  values  of  relaxation 
result  in  smaller  value  differences  per  iteration.  However, 
excessively  small  relaxation  values  require  additional 
computer  time.  A  balance  gained  by  user  experience  enables 
optimal   convergence  with  minimal  computer  time. 

When  multiple  runs  are  required  for  large  programs  output 
data  is  stored  in  a  "DF09"  file  which  allows  the  next  run  to 
resume  the  computations   (i.e.,  a  restart)  where  they  last 
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ended.  The  coiriplexity  of  the  problem  and  the  required 
accuracy  determine  the  quantity  of  computer  runs  needed.  The 
number  of  runs  can  vary  from  the  tens  to  the  thousands  and 
higher . 

B-     ASSUMPTIC»>JS     HSr    THE    MODEX. 

There  are  several  approximations  that  deserve  mention  in 
the  modeling  of  the  Naval  Air  Station,  Jacksonville,  Florida 
Hush  House.  Because  the  geometry  of  the  model  uses  cartesian 
coordinates,  the  actual  Hush  House's  rounded  and  curved 
surfaces,  including  the  augmenter  tube,  Hush  House  overhead, 
engine  and  exhaust  gas  exit  ramp,  were  approximated  as  flat 
stepped  surfaces. 

The  inlet  "muffler"  arrangement  has  been  eliminated  for 
simplification.  The  corrugated  metal  found  throughout  the 
structure  to  cover  the  sound  absorption  material  was 
considered  to  have  a  friction  surface.  The  model,  therefore, 
included  friction  for  all  walls,  floors,  etc. 

The  actual  test  data  used  for  comparison  between  the 
model  and  facility  contain  several  notable  differences.  The 
model  uses  the  symmetry  of  the  centerline  of  the  Hush  House 
to  reduce  the  number  of  cells  by  half.  The  engine  is  located 
on  the  centerline  in  an  ideal  cradle  position.  Because  of 
the  symmetry  only  half  of  the  engine  is  observed  as  if  the 
Hush  House  was  being  used  as  a  test  cell  without  the  aircraft 
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present.  Since  the  correct  engine  exhaust  area  is  used  the 
units  of  mass  flow  rate  and  temperature  are  correct  for  the 
model  configuration.  The  engine  is  configured  for  afterburner 
operation.  The  physical  aircraft  is  not  present. 

Actual  test  data  used  for  comparison  consists  of  the 
entire  aircraft  being  placed  in  the  facility.  The  engine 
configuration  has  both  engines  running  with  one  at 
afterburner  while  the  other  engine  is  at  "engine  idle".  The 
obvious  result  of  the  engine  arrangement  is  a  non-symmetrical 
engine  location  in  comparison  to  the  symmetric  model. 

c-    mode^jItjg   the   PR0BUE3^ 

Data,  including  blueprints  of  the  Hush  House  and  aircraft 
engine  information,  was  gathered  to  accurately  depict  the 
facility  configuration.  The  information  obtained  was  in 
English  units  and  was  converted  to  the  International  Metric 
System  or  "S.I."  using  standard  conversion  factors. 

With  the  appropriate  information  gathered  it  was  time  to 
"fit"  the  Hush  House  to  the  cartesian  coordinate  system  of  X, 
Y  and  Z  coordinates  using  PHOENICS  version  1.3.  The  cartesian 
coordinate  system  was  choosen  because  of  its  simplicity.  Since 
the  actual  structure  does  have  rounded  surfaces  some 
approximations  were  made  specifically  to  the  augmenter  tube 
and  overhead  of  the  inlet  bay.  Additionally,  the  aircraft  was 
eliminated  to  reduce  the  complexity  of  the  problem  allowing 
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the  structure  to  serve  in  its  "Test  Cell"  configuration. 
Drawings  are  found  in  Figures  3.2  and  3.3  of  this  section 
illustrating  the  key  items  found  in  the  structure.  The 
drawings  view  through  the  structure  to  note  the  key  items, 
such  as:  air  inlets,  augmenter  tube  and  engine  locations. 
Additionally,  the  structure  walls  were  considered  to  be  of 
infinitesimal  thickness  to  further  simplify  the  grid 
structure.  The  engine  was  simplified  from  a  cylindrical  shape 
to  one  of  a  rectangle  of  equal  intake  and  exhaust  area  as  the 
original  cylinder.  Other  small  items,  such  as  the  engine 
stand,  were  eliminated. 

To  take  advantage  of  the  symmetry  of  the  Hush  House  and 
thereby  reduce  computer  CPU  time,  only  half  of  the  Hush  House, 
engine  and  augmenter  tube  were  modelled  in  the  "Ql"  file.  This 
was  easily  accomplished  by  dividing  the  structure  at  the 
centerline.  This  also  appreciably  reduces  the  size  of  the 
"Ql"  file. 

Figures  3.2  and  3.3  graphically  display  the  X-Z  and  the 
Y-Z  planes  of  the  structure.  Intuition,  trial-and-error ,  and 
available  computer  time  dictated  the  grid  spacing  used. 
Obviously,  volume  areas  of  interest  such  as  the  engine  exhaust 
and  augmenter  required  smaller  cell  size. 

The  "Ql"  file  describes  the  structure  using  several  types 
of  statements.  The  "CONPOR"  statement  allows  the  user  to 
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define  walls,  overheads  and  the  floor.  The  "PATCH"  and 
"COVAL"  statements  work  together  to  indicate  items  such  as 
an  engine  exhaust  or  atmospheric  boundaries.  Dependent 
variable  values  for  say  the  atmospheric  boundaries  can  be 
"fixed"  or  allowed  to  "float"  depending  on  the  problem 
definition.  When  the  walls,  inlet  bay  overhead,  augmenter 
tube,  engine  and  floor  dimensions  have  been  measured,  a  grid 
for  the  system  can  be  constructed,  based  on  locations  and 
anticipated  responses  over  different  "areas"  of  the  system. 
Areas  of  anticipated  large  response  gradients  would  have 
smaller  cell  volumes  than  low  activity  areas. 

Note  that  in  figure  3.2  the  Z-axis  consists  of  40  cells 
consecutively  numbered  from  1  to  40  going  left  to  right.  The 
X-axis  consists  of  10  cells  consecutively  numbered  1  through 
10  going  from  lower  tc  upper.  Figure  3.3  is  similarly  numbered 
where  the  Y-axis  runs  frorr.  1  to  12  going  from  lower  to  upper 
and  the  Z-a:-:is  is  the  same  as  previously  described. 

In  addition  to  describing  the  geometric  structure  and 
prescribing  the  "RELAXATION"  devices,  the  properties  of  the 
fluid  medium  must  be  specified.  The  following  boundary 
conditions  were  used: 

1.  ambient  temperature  =  295  K 

2.  ambient  pressure  =  101325  Pa 

3.  gas  constant  =  286  Joules/ (kg  -  K) 
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4.  Cp  =  1004  Watt-sec/(kg  -  K) 

5.  mass  flow  rate  (engine)  =  197.72  kg/(sec-sq  meters) 

6.  temperature  (engine  out)  =  1944  K 

Additionally,  "friction"  was  applied  to  all  surfaces  of 
the  Hush  House.  Appendix  A  contains  the  final  form  of  the 
Ql  file  for  reference.  An  in-depth  explanation  of  all  the  "Ql" 
program  parameters  can  be  found  in  the  "Beginners  Guide  and 
Reference  Manual"  (See  Reference  4). 

Several  methods  and  theories  have  been  previously  used  to 
attempt  to  reach  a  "converged"  solution  in  other  test 
facilities.  The  method  chosen  for  this  project  was  to 
initially  minimize  the  "excitation"  of  the  model  thereby 
reducing  the  "jump"  required  in  solving  the  momentum, 
enthalpy,  turbulent  kinetic  energy  dissipation  rate,  turbulent 
kinetic  energy  and  mass  equations.  To  accomplish  this,  during 
successive  runs  (where  a  "run"  consists  of  several  "sweeps"), 
small  incremental  increases  to  the  source  were  applied,  i.e., 
jet  engine  exhaust  temperature  and  mass  flow  rate  were  slowly 
increased,  frorr;  run  to  run,  while  carefully  monitoring  the 
residuals  as  well  as  the  mass  flow  rate  in  and  out  of  the 
structure.  The  turbulence  model,  "KEMODL",  was  used 
throughout  rather  than  considering  the  flow  first  as  laminar 
and  then  adding  turbulence.  This  method  appears  to  have  met 
with  considerable  success.  The  table  that  follows  illustrates 
the  iterative  process  used  throughout  the  computer  modeling: 


Table  3.1   ITERATIVE  PROCESS 


Temperature 


Flow  Rate 


Sweeps 


400 
500 
750 

750 
1000 
1250 
1250 
1250 
1500 
1750 
1750 
1944 
1944 
1944 


10 
25 
25 

50 

75 

75 
100 
125 
125 
125 
150 
150 
175 
197.72 


600 
600 
600 

600 
2500 
1200 
1200 
1200 
1200 
2000 
2000 
2000 
5000 
5500 


Total  Sweeps  =  26,200 
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This  amounts  to  approximately  147  hours  of  computer  CPU  time 
on  the  IBM  3033  mainframe  computer  installed  at  NPS. 

Using  the  MVS  "batch"  system  on  the  IBM  3033  mainframe 
proved  to  be  time  consuming.  The  PHOENICS  Code  is  designed  to 
work  both  interactively  and  in  batch  mode.  Due  to  the  high 
computer  usage  at  NPS  the  "interactive"  mode  was  not 
available.  Limited  to  the  batch  mode,  careful  planning  for 
evening  or  night  runs  was  required  to  maximize  CPU  time.  The 
number  of  "Degrees  of  Freedom"  (DOF)  that  must  be  solved  for 
gives  an  idea  of  the  problem  complexity.  DOF  =  NX  *  NY  *  NZ 
*  the  7  dependent  variables.  Solving  this  for  the  grid  size 
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of  10  *  12  *  41  yields  nearly  35,000  values  calculated  per 
each  sweep.  Typically,  150  to  200  sweeps  could  be  accomplished 
during  the  day  using  four  15  minute  (class  "G")  runs  of  50 
sweeps  each.  After  normal  working  hours  three  one  hour  (class 
"J")  runs  totalling  500  to  600  sweeps  were  accomplished.  Data 
stored  on  the  all  important  DF09  file  was  continuously 
"backed-up"  after  a  successful  run  to  prevent  loss  of  data  if 
an  improper  subsequent  run  resulted  in  erroneous  data.  The 
backed-up  data  could  then  be  restored  to  the  DF09  file  to 
enable  continuation  of  the  runs  without  starting  at  the  very 
beginning . 
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Figure  3.2   Hush  House  Model  (X-Z  plane) 
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Figure  3.3   Hush  House  Model  (Y-Z  plane) 


26 


XV  -     RESULTS 

A-     THE    FINAL    SOLLmOISr 

After  obtaining  the  "solution"  for  the  problem  the 
question  still  remaining  is  just  how  "good"  are  the  computer 
model  results?  Some  of  the  "tools"  available  for  use,  such  as 
checks  on  residual  values  and  mass  flow  rate  calculations, 
were  discussed  in  Chapter  III.  The  checks  indicate  that  a 
converged  solution  is  obtained,  but,  still  do  not  prove  the 
computer  model  data  validity  with  the  actual  facility  data. 

Using  actual  (experimental  field  measurements)  results 
from  Reference  2  allows  for  an  accurate  evaluation  of  the 
computer  model.  Appendices  "C"  and  "D"  offer  some  experimental 
data  for  velocities  and  temperatures,  thereby  indicating  mass 
flow  rate  and  augmentation  ratio.  Additionally,  figures  4.1, 
4.2  and  4.3  found  in  this  section  give  the  numerical 
(computer)  model  results  of  differential  pressure,  temperature 
and  axial  velocity  profiles  encountered  at  the  IX=1  slab. 
Graphics  generated  pictures  (generated  by  "PHOTON"  which 
converts  the  computer  generated  results  to  graphical  format) 
allow  for  the  interpretation  and  digestion  of  the  nearly  300 
pages  of  numerical  output  of  a  PHOENICS  run.  A  condensed 
output  is  found  in  Appendix  "B".  In  this  case,  PHOTON  graphics 
specifically  provides  graphical   information  for  a  better 
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interpretation  of  velocity  fields  and  both  pressure  and 
temperature  profiles  and  contours  enabling  a  better 
understanding  of  computer  generated  data.  Use  of  all  these 
additional  "tools"  allow  for  a  determination  as  to  whether  the 
computer  results  are  valid. 

Observing  figures  4.4  and  4.5  (top  looking  down)  show  the 
flow  pattern  of  the  air  entering  the  five  air  inlets  into  the 
Hush  House  at  IY=4  (engine  level).  A  "typical"  low  velocity 
air  flow  region  (Figure  4.4),  illustrated  by  the  minuscule 
length  of  the  velocity  arrows  (velocities  of  less  than  2 
meters/sec  or  5.5  feet/sec),  is  noted  between  the  first  two 
and  final  three  inlets.  Figures  4.6,  4.7  and  4.8  illustrate 
the  side  view  velocity  flow  patterns  for  the  overall  facility, 
engine  area,  and  exhaust  to  atmosphere  respectively.  Note  the 
ever  reducing  velocity  magnitudes  which  vary  from  818  meters 
per  second  (2684  feet/sec)  at  the  engine  exhaust  to  less  than 
100  meters  per  second  (328  feet/sec)  as  the  gases  travel 
further  away  from  the  engine  and  down  the  augmenter  tube  until 
reaching  the  exhaust. 

Figures  4.4,  4.5  and  4.9  (top  looking  down)  show  the  air 
and  gas  flow  patterns  entering  and  exiting  the  engine  region 
and  augmenter  tube  at  engine  level.  Figure  4.10  indicates  the 
flow  patterns  at  the  top  of  the  augmenter  tube  (IY=8)  for  the 
entire  Hush  House.  Of  interest  is  the  counter-flowing  currents 
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at  the  top  of  the  augmenter  tube  starting  at  the  engine 
exhaust  and  running  to  nearly  half  way  down  the  augmenter  tube 
itself . 

The  engine  occupies  cells  defined  by  IX=1 ,  IY=4  to  5  and 
IZ=12  to  14.  As  anticipated,  figures  4.4,  4.9  and  4.11  show 
that  at  engine  level  the  air  in  the  Hush  House  moves  toward 
the  engine  inlet  from  the  surrounding  area  while  air  also  is 
entrained  at  the  exhaust  of  the  engine.  Figure  4.1  illustrates 
this  entrained  "effect"  noted  by  the  high  differential 
pressure  of  the  exhaust  gases  immediately  exiting  the  engine 
at  IX=1  and  IY=4.  Less  than  a  meter  after  exiting  the  engine 
exhaust,  the  exhaust  gases  and  inlet  air  have  begun  mixing 
which  results  in  an  area  of  drastic  temperature  reduction  for 
the  exhaust  gases  in  the  Hush  House  and  immediate  augmenter 
area.  Gradually,  as  the  gases  pass  further  downstream, 
pressure  begins  rising  which  reduces  the  velocity  of  the  gases 
in  the  augmenter  tube  until  the  gases  exit  to  the  atmosphere. 
At  the  inlet  area  to  the  engine  (at  Z=18  meters  on  figure  4.1) 
a  low  pressure  area  is  noted  as  one  would  also  anticipate  as 
the  air  is  drawn  into  the  engine  inlet. 

The  exhaust  gases  exit  the  engine  at  1944  Kelvin  (3040  F) . 
The  inlet  air  and  exhaust  gases  mix  due  to  the  entrainment 
"effect"  causing  a  distinct  drop  (See  figure  4.2)  in 
temperature  as  the  gases  travel  further  down  the  augmenter 
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tube.  Temperature  contours  at  the  IX=1  slab  shown  in  figures 
4.12  and  4.13  portray  the  temperature  patterns  out  of  the 
engine  progressing  down  and  out  of  the  augmenter  tube. 
Analogous  temperature  contours  at  the  IX=3  slab  are  presented 
in  figures  4.14  and  4.15.  Figures  4.16,  4.17  and  4.18  embody 
the  temperature  contours  at  the  IY=5  (engine)  level  while 
figure  4.19  shows  the  temperatures  at  the  top  of  the  augmenter 
tube  at  IY=8.  The  temperature  contours  indicate  a 
progressively  widening  area  of  decreasing  exhaust  gas 
temperatures  where  the  vast  majority  of  the  temperature  drop 
occurs  in  the  first  five  meters  after  the  engine  exhaust  (the 
temperature  differential  goes  from  1944  K  to  800  K) . 

Differential  pressure  contours  of  the  engine  exhaust  at 
IX=1  are  shown  in  figures  4.20  (IZ=  14  to  24)  and  4.21  (IZ= 
25  to  40) .  As  indicated  in  the  pressure  profile,  the  exhaust 
jumps  to  nearly  5000  pascals  above  ambient  upon  exiting  the 
engine,  but  soon  after  drops  into  a  low  pressure  region  until 
reaching  ambient  pressure  half  way  down  the  augmenter  tube  at 
IZ=24.  Gradually,  the  pressure  builds  until  the  gases  exhaust 
to  the  atmosphere.  Similar  contours  for  the  IX=3  slab  are 
included  in  figures  4.22  and  4.23.  Pressure  contours  for  the 
X-Z  planes  at  IY=5  (engine  level)  and  IY=8  (top  of  the 
augmenter  tube)  are  contained  in  figures  4.24,  4.25  and  4.26. 
Of  noted  interest  is  how  the  differential  pressures  remain 
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fairly  uniform  at  a  given  IZ  location  as  the  gases  travel  down 
the  augmenter  tube  with  the  exception  of  the  region  initially 
downstream  from  the  engine  exhaust. 

B-  ACTUAL  VERSUS  MODEHL.  ODMRARISOMS 

Generally,  the  results  from  the  computer  model  closely 
approximated  the  actual.  The  below  "actual"  data,  and  that  of 
Appendix  "C"  and  "D"  for  the  computer  model  are  indicated  in 

both  "English"  and  "S.I."  units  below.  A  comparison  of  the 
PHOENICS  output  visually  shown  in  Appendices  "C"  and  "D" 
illustrates  the  dissimilarity  between  the  actual  and  computer 
model.  As  indicated  by  the  data,  the  differences  in  the 
computer  model  (235.5  ft/sec  or  71.8  meters/sec)  and  actual 
(227.2  ft/sec  or  59,3  meters/sec)  averaged  velocities  found 
in  the  augmenter  tube  at  slab  IZ=27  differed  by  less  than  4 
percent.  Averaged  temperatures  in  the  same  location  differed 
between  401.6  degrees  Fahrenheit  (478  K)  for  the  computer 
model  and  454  degrees  Fahrenheit  (507  K)  actual  with  a 
resultant  difference  of  6  percent.  Reference  2  does  not 
indicate  whether  the  effects  of  radiant  energy  from  the  engine 
exhaust  are  considered  for  the  actual  temperature  readings 
taken  in  the  augmenter  tube.  If  these  effects  were  not 
considered  the  6  percent  difference  would  be  reduced.  The  Hush 
House  mass  flow  rate  for  the  computer  model  (2373  Ibm/sec  or 
1076  kg/sec)  versus  the  actual  (2067  Ibm/sec  or  938  kg/sec) 
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differed  by  14.8  percent.  The  computer  model  mass  flow  rate 
value  was  actually  doubled  (hence  the  2373  Ibm/sec  or  1076 
kg/sec  value)  to  "model"  the  entire  engine  and  Hush  House. 
The  Augmentation  Ratio  is  defined  as  follows: 


Mlnlet   -  Menglne 

A.R.  =  

Mlnlet 


where : 

Mlnlet  =  inlet  mass  flow  rate 
Mengine  =  engine  mass  flow  rate 

The  inlet  mass  flow  rate  =  air  density  *  five  inlet 
opening  areas  *  averaged  inlet  velocity.  The  engine  mass  flow 
rate  is  given  at  197.72  kg/sec/sq  meters  *  0.195  sq  m.eters 
engine  exhaust  area  for  half  of  the  engine. 

Augmentation  ratios  range  from  5  to  15  depending  on  the 
design  of  the  Hush  House  or  Test  Cell  where  the  higher  the 
ratio  typically  indicates  a  more  efficient  design.  The 
computer  model  augm.entation  ratio  was  found  to  be  10.6  while 
the  actual  was  8.3.  The  resultant  difference  is  therefore 
nearly  28  percent.  Considering  the  limited  data  taken  at  this 
actual  power  level  (See  Appendix  C)  it  is  quite  conceivable 
that  the  actual  augmentation  ratio  could  be  greater. 

For  the  actual  results  the  velocity,  pressure  and 
temperature  measuring  devices  were  photographed  because  of  the 
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large  fluctuations  of  the  gauges.  Additionally,  velocity  data 
at  the  final  three  inlets  could  not  accurately  be  taken  making 
it  difficult  to  accurately  determine  the  augmentation  ratio 
and  the  true  velocities  at  the  inlets. 
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Figure  4.1 


Differential  Pressure  Profile  at  Volume  Slab 
IX=1  for  IY=4  (Engine  Height) 
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Figure  4.2   Temperature  Profile  at  Volume  Slab  IX=1 
for  IY=4  (Engine  Height) 
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Figure  4.3   Axial  Velocity  Profile  at  Volume  Slab  IX=1 

for  IY=4  (Engine  Height) 
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Figure  4.4   Flow  Vector  Plot  of  Hush  House  Inlets  and 
Engine  Area  (X-Z  Plane  at  IY=4) 
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Figure  4.5   Flow  Vector  Plot  in  X-Z  Plane  at  IY=4 
(Engine  Exhaust  and  Augmenter  Tube) 
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Figure    4.6 


Flow  Vector  Plot  in  Y-Z  Plane  at  IX=1 
(Entire  Grid) 


39 


® 


>  <- 


N 


X 


LU 

Z 

a 

_j 

Q. 

N 
I 

> 


CD 
IX 
O 

t— 

u 

LU 
> 


u 
o 
_1 
IJJ 
> 


L 
O 

u 

(D 


E 
O 
O 

r\i 


A 


A 


u 

z 

UJ 

o 

X 
CL 


LU 

m 

z 
a 

I 
u 

I 
Ln 

Z) 

I 

LU 


> 

Z 

o 
tn 
i<: 
u 
a 
•-) 

LD 

cr 


Figure    4.7 


Flow  Vector  Plot  in  Y-Z  Plane  at  IX=1 
(Engine  Exhaust  Area) 
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Figure    4 . 8 


Flow  Vector  Plot  in  Y-Z  Plane  at  IX=1 
(Exhaust  to  Atmosphere) 
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Figure  4.9 


Flow  Vector  Plot  in  X-Z  Plane  at  IY=4 
(End  of  Augmenter  Tube) 
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Figure  4.10   Flow  Vector  Plot  in  X-Z  Plane  at  IY=8 
(Mirrored  Image  of  Entire  Grid) 
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Figure  4.11 


Flow  Vector  Plot  of  Augmenter  Tube  Area 
for  X-Z  Plane  at  IY=4 
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Figure  4.12   Temperature  Contour  in  Y-Z  Plane  at  IX=1 
(Engine  and  Augmenter  Tube) 
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Figure  4.13   Temperature  Contour  in  Y-Z  Plane  at  IX=1 
(Augmenter  Tube  and  Exit) 
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Figure    4.14      Temperature    Contour    in   Y-Z    Plane    at    IX=3 

(Engine    Exhaust) 
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Figure  4.15   Temperature  Contour  in  Y-Z  Plane  at  IX=3 
(Engine  Exhaust  to  Atmosphere) 
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Figure  4.16   Temperature  Contour  in  X-Z  Plane  at  IY=5 
(Mirrored  Image  of  Engine  Exhaust) 
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Figure  4.17   Temperature  Contour  in  X-Z  Plane  at  IY=5 

(Engine  Exhaust) 
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Figure  4.18   Temperature  Contour  in  X-Z  Plane  at  IY=5 
(Engine  Exhaust  in  Augmenter  Tube) 
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Figure  4.19   Temperature  Contour  in  X-Z  Plane  at  IY=8 
(Engine  Exhaust  in  Augmenter  Tube) 
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Figure    4.20 


Differential  Pressure  Contour  in  Y-Z  Plane 
at  IX=1  (Engine  Exhaust) 
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Figure  4.21   Differential  Pressure  Contour  in  Y-Z  Plane 
at  IX=1  (Engine  Exhaust  to  Atmosphere) 
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Figure    4.23 


Differential  Pressure  Contour  in  Y-Z  Plane 
at  IX=3  (Exhaust  to  Atmosphere) 
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Figure    4.24 


Differential  Pressure  Contour  in  X-Z  Plane 
at  IY=5  (Engine  Exhaust) 
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Figure  4.25   Differential  Pressure  Contour  in  X-Z  Plane 
at  IY=5  (Augmenter  Tube) 
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Figure  4.26   Differential  Pressure  Contour  in  X-Z  Plane 

at  IY=8  (Augmenter  Tube) 
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V-      SUMMARY    AND    RECOMMEaSTDATICSSrS 

A-     SUMMARY 

Using  the  Parabolic,  Hyperbolic,  or  Elliptic  Numerical 
Integration  Code  Series  (PHOENICS)  a  computer  model  of  the 
Naval  Air  Station  Jacksonville  Florida's  Hush  House  facility 
with  a  J-79-GE-8  gas  turbine  engine  was  simulated.  The 
performance  of  this  computer  model  was  compared  to  actual 
performance  data  taken  with  an  F-4  Phantom  II  undergoing 
testing  at  an  identical  facility  located  at  Kelly  AFB .  An  on- 
site  inspection  was  conducted  at  Dobbins  AFB  in  Atlanta 
Georgia  to  provide  "hands-on"  experience  with  the  facility. 

The  computer  analysis  consisted  of  slowly  building  up 
excitation  in  the  model  (using  the  KEMODL  turbulence  model) 
until  full  power  was  obtained.  Computer  generated  (tabulated 
and  graphical)  results  were  compared  with  the  actual  test 
data.  Results  of  the  computer  model,  in  comparison  to  the 
actual  data,  were  in  good  agreement.  Flow  directions, 
velocities  and  temperatures  were  consistent  with  actual 
results . 

B .  RECX^lMEasrDATIOiSrS 

A  comparison  of  actual  to  computer  model  data  for  this 
thesis   indicates   that   PHOENICS   is   a   useful   design   and 
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troubleshooting  tool  allowing  the  user  to  better  understand 
the  physics  of  the  model.  PHOENICS  at  NPS  is  now  just  coining 
into  potent  use  with  a  variety  of  research  projects  for 
aeronautical  and  mechanical  engineers. 

To  ensure  proper  and  full  use  of  this  code  several  major 
"glitches"  must  be  overcome.  Time  constraints  placed  on  the 
use  of  the  IBM  3033  main-frame  computer  limit  the  number  of 
"sweeps"  obtainable  by  the  student.  The  batch  system  does  not 
lend  itself  to  correcting  possible  errors  that  may  arise.  An 
interactive  system,  such  as  the  "VAX"  stand  alone  systems 
being  installed  in  Mechanical  Engineering,  is  an  ideal  answer. 

PHOTON  was  not  available  with  the  PHOENICS  installation 
on  the  IBM  3033  main-frame  computer.  A  complicated  conversion 
process  (See  Appendix  "E")  was  undertaken  to  convert  IBM  to 
VAX  "format"  to  allow  the  use  of  PHOTON  using  data  from  the 
main-frame.  This  conversion  process  which  involved  a  swap  of 
information  from  IBM  to  DEC  to  VAX  hardware  took  nearly  2 
weeks  to  accomplish.  In  September  1988  this  transfer  procedure 
will  no  longer  be  possible  at  NPS  with  the  elimination  of 
certain  hardware  at  the  Computer  Science  Department  although 
data  transfer  will  still  be  possible  to  other  facilities  such 
as  CHAM. 

The  use  of  PHOTON  early  in  the  computer  "runs"  would  have 
uncovered  several  missing  "CONPOR"  statements  and  therefore 
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"verified"  the  structural  integrity  of  the  facility. 
Additionally,  the  amount  of  time  required  to  "wade"  through 
the  300  pages  of  printout  could  be  reduced  to  several  video 
screens  with  a  much  reduced  volume  of  printout. 

The  computer  model  used  in  this  thesis  was  formulated 
using  cartesian  coordinates.  Several  parts  of  the  structure, 
including  the  augmenter  tube  and  engine,  were  curve  shaped. 
Geometric  approximations  were  made  to  substitute  for  the  curve 
shaped  areas.  A  more  computationally  time  consuming  model  with 
increased  numbers  of  cells  and/or  the  use  of  body-fitted- 
coordinates  (BFC)  would  enhance  the  computer  model. 

As  noted  in  Figure  4.2,  over  fifty  percent  of  the 
reduction  in  engine  exhaust  gas  temperatures  occurred  in  the 
first  20  percent  of  the  augmenter  tube.  Could  an  improved 
design  include  shortening  the  augmenter  tube?  This  type  of 
"what  if"  question  can  be  answered  using  PHOENICS .  With  the 
variable  data  file  ("DF09")  already  in  place  the  grid  pattern 
can  be  modified  to  shorten  the  augmenter  tube.  Minimal 
computations  would  be  required  to  gather  definitive  results. 

Parameter  studies,  including  shortening  the  length  of  the 
augmenter  tube  (to  determine  the  effect  of  augmenter  tube 
length  on  performance)  ,  movement  of  the  engine  location  in  the 
Hush  House,  changing  the  inlet  locations  and  sizes,  and  the 
use  of  different  turbulence  models  could  lead  to  improvements 
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on  the  present  Hush  House  design.  In  fact,  the  ability  to 
change  parameters  and  observe  their  effect  on  systein 
performance  is  the  prime  attribute  of  numerical  modeling. 

Finally,  a  discussion  of  turbulence  is  in  order.  PHOENICS 
offers  four  turbulence  models  which  can  be  readily  selected. 
The  four  models  are:  Algebraic,  Mixing  Length,  K-Mixing  Length 
or  K-Epsilon.  The  user  can  also  formulate  a  model  of  his  own 
choosing  which  can  be  incorporated  into  the  code  via  the 
"GROUND"  program.  Two  questions  arise.  First,  which  turbulence 
model  is  appropriate  for  the  problem  under  investigation? 
Secondly,  how  strong  is  the  coupling  between  other  system 
response  variables  and  turbulence  variables?  The  second 
question  might  be  answered  by  running  the  PHOENICS  Code  on  the 
problem  using  several  different  turbulence  models. 
Unfortunately  there  was  not  adequate  time  to  perform  this  in 
the  investigation.  The  first  question  can  not  be  resolved 
until  a  better  understanding  of  turbulence  is  obtained. 

C .  CCfNCUUSXCfNS 

The  computer  model  results  generated  by  this  thesis 
closely  parallel  the  actual  test  data.  The  results  therefore 
indicate  that  the  PHOENICS  code  can  be  successfully  utilized 
to  model  future  Test  Cell  and  Hush  House  applications  as  well 
as  "trouble-shoot"  suspect  areas  of  interest.  Also,  areas  that 
may  be  difficult  to  properly  monitor  with  test  equipment  could 
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be  modelled  instead.  With  the  use  of  additional  "cells"  and/or 
BFC  a  more  refined  model  could  more  closely  duplicate  the  test 
facility  providing  more  accurate  results.  Follow-up  parametric 
studies  to  improve  and  verify  the  design  and  data  seem  in 
order.  The  results  indicate  that  the  Hush  House  is  performing 
as  intended. 
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APPE3MDIX    A 
THE    PHOENXCS    Ql    FILE 


TALK=F,RUN(1,1) 

•  NAVAL  AIR  STATION  JACKSONVILLE  FLORIDA  HUSH  HOUSE  TEST  « 
«  FACILITY.  THESIS  PROJECT  OF  LT  ERIC  A.  NICOLAUS  USCG    « 

•  WITH  THESIS  INSTRUCTOR  PROF  DAVID  SALINAS.  • 

W)rXK»K>()(»>>)(it¥l«:iMW)C)CMK«#>CMmii(m<>l)mit<<iCltlli<MKKKXK)>M)(ltllMI<M)(l>ICIt« 

»                LIST  OF  VARIABLES  » 

«  « 

»  HJET:  VELOCITY  OF  ENGINE  EXHAUST  JMETERS/SEC)  » 

»  MJET:  MASS  FLOH  RATE  OF  GAS  ( KG/SECOND/SQ  METER)  » 

«  TJET:  JET  EXHAUST  GAS  TEMPERATURE  (DEGREES  KELVIN)  « 
«  PJET:  EXHAUST  PRESSURE  OF  JET  ENGINE  (NEHTONS/SQ  METER)* 

«  RHOJET:  DENSITY  OF  THE  GAS  (KG/CUBIC  METER)  » 

«  HJET:  ENTHALPY  (JOULES/KG  OR  NEWTON»METERS/KG )  » 

«  KEJET:  KINETIC  ENERGY  (JOULES  OR  NEHTON-METERS )  « 

«  GASCON:  GAS  CONSTANT  ( JOULES/(  KG-DEG  KELVIN))  « 

«  TAMB:  AMBIENT  AIR  TEMPERATURE  (DEGREES  KELVIN)  « 

«  HAKB:  AMBIENT  ENTHALPY  (JOULES/KG)  » 

«  KEINJ:  TURBULENT  KINETIC  ENERGY  (JOULES)  OF  JET  « 

«  CSUBP:  SP  HT  3  CONST  PRE  ( JOULES/C KG-DEG  KELVIN))  « 

•  EPINJ:  DISSIPATION  RATE  OF  TURBULENT  KE  (JOULES)  OF  JET« 
»  KEINA:  TURBULENT  KINETIC  ENERGY  (JOULES)  AMBIENT  » 
»  EPINA:  DISSIPATION  RATE  OF  TURBULENT  KE  (JOULES)  AMB  « 
«  PI:  PRESSURE  (NEWTONS/SCJUARE  METER)  • 
«  HI:  ENTHALPY  (JOULES/KG)  « 

•  Ul:  VELOCITY  IN  X-DIRECTION  (METERS/SEC)  « 
«  VI:  VELOCITY  IN  Y-DIRECTION  (METERS/SEC)  » 
«  Wl:  VELOCITY  IN  Z-DIRECTION  (METERS/SEC)  « 
«  TMPl:  TEMP  OF  THE  FIRST  PHASE  (DEGREES  KELVIN)  » 

•  RHOl:  DENSITY  OF  THE  FIRST  PHASE  (KG/CUBIC  METER  I  • 
«  PRESSO:  REFERENCE  ATMOSPHERIC  PRESSURE  ( NEWTOtJS/StJ  mtr)« 
»  TMPIA:  GROUrJD  REF  TEMP  USED  IN  GRND2  (  DEG  KELVIN)  • 
«  TMPIB:  1/CSUBP  USED  IN  GRM02  EC3U  (KG-DEG  KE LVIN/ JOULES )» 
«  DRHIDP:  DELTA  RHO/DELTA  PRESSURE  H/RESPECT  TO  HI  • 
«  FACT:  CORRECTION  OR  "FUDGE"  FACTOR  « 
«  ENUL:  LAMINAR  KINEMATIC  VISCOSITY  (METERS  SQ/SEC )  » 
«  ENUT:  TURBULENT  KINEMATIC  VISCOSITY  (METERS  SQ/SEC)  « 

GROUP  1.  RUN  TITLE  AND  OTHER  PRELIMINARIES 
TEXTINAS  JACKSONVILLE  HUSH  CHAMBER) 

REAL(  HJET, MJET, TJET, PJET, RHOJET, HJET, KEJET, GASCON, CSUBP, TAMB) 
REAL(HAMB, KEINJ, EPINJ, KEINA, EPINA) 

GROUP  2.  TRANSIENCE)  TIME -STEP  SPECIFICATION  ( DEFAULT=STEADY ) 

GROUP  2.  X-DIRECTION  GRID  SPECIFICATION 
NX  =  10 

XFRACI1)=-1  }  XFRAC(2)=0.21 
XFRAC(3)='+   i  XFRACI^)  =  0.5'+7S 
XFRAC(5)=3   )  XFRAC16)=2.0 
XFRAC17)  =  2   )  XFRACl  8  )  =  !.<+ 
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GROUP  4.  Y-DIRECTION  GRID  SPECIFICATICN 


NY=12 

YFRAC(1)=-S 
YFRACI5)=2 
YFRACI S)=2 
YFRAC(7)=<f 

GROUP  5. 

2FRAC(1)=-1 

2FRAC(3)=1 

2FRACt5)=l 

2FRAC(7J=1 

ZFRAC(9)=2 

ZFRAC(11)=1 

2FRAC113)=1 

2FRAC(1S)=1 

ZFRAC(a7)=2 

2FRACJ19)=1 

2FRAC(21)=2 

2FRAC(23)=5 

2FRAC(25)=2 

2FRAC(27)=7 

2FRAC(29]=1 

2FRAC(31)=^ 

ZFRAC(33)=<+ 

ZFRAC(3S)='i 

2FRAC(37)=1 


i  YFRAC(2J=0.S62 

>  YFRAC(4)=0.31't 

>  YFRAC(6)=0.562 
\  YFRAC(8)=0.7S 

2-DIRECTION  GRID  SPECIFICATION 


ZFRACI 
2FRAC( 
2FRAC1 
2FRAC( 
2FRAC( 
2FRAC( 
2FRAC( 
2FRACI 
2FRAC( 
2FRAC( 
2FRAC( 
2FRAC( 
2FRAC( 
2FRAC( 
ZFRACt 
2FRAC( 
ZFRAC( 
ZFRACC 
ZFRAC( 


2)=1.0 

^)=1.7 

6)  =  2.0 

8)=1.7 

10)=3.0 

12)=1.S 

14)=1.0 

16)=1,7 

18)=1.0 

20)=1.7 

22)=1.0 

2<*)-0.S66666 

2fe)=0.3S 

28)=3.0 

30)=2.0 

32)=0.7S 

3^)=0.S7S 

36  )=0.5 

38  J =3.0 


GROUP  6.  BODY-FITTED  COORDINATES  OR  GRID  DISTORTION 

GROUP  7.  VARIABLES  STORED,  SOLVED  &  NAMED 
S0LVE(P1,H1,U1,V1,H1) 
S0LUTNIP1,Y,Y,Y,N,N,N) 
ST0REITMP1,RH01) 

GROUP  8.  TERMS  (IN  DIFFERENTIAL  EQUATIONS)  t   DEVICES 

GROUP  9.  PROPERTIES  OF  THE  MEDIUM 
TAMB=295. 
TJET  =  19'+'+. 
PJET=0. 
MJET=197.72 
PRESS0=101325. 
RH01A=0. 
GASC0N=286. 
RHO IB =1. /GASCON 

RHOJET=((PRESSO*PJET)«RH01B)/TJET 
HJET=MJET/RHOJET 
TMP1A=TINY 

»*«  TMPIB  =  1.  DIVIDED  BY  C-SUB-P)  USED  IN  TMP1=GRND5 
CSU3P  =  100'+. 

hamb=tamb»csubp 

tmpie;=i./csubp 

hjet=csubp«tjet 
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»*   NEXT  2  STATEMENTS  INVOKE  COMPRESSIBLE  FLOW  IN  CREXl 
TMP1=GRND2 
RH01=GRNDS 
DRH10P-6RN05 
ENUL=l.E-05 
ENUT=200.»ENUL 
TURMOOIKEMODLJ 
STORE (VIST) 

KEINJ=0.5*(0.00S»WJET)**2  >  EPINJ=0.09»KEINJ*«l.S/( 0.05*0.5  ) 
KEINA=0.5«l .005*1.0 )«*2  )  EPINA=0.09«KEINA««l.S/( 0. 1 ) 

GROUP  10.  INTER-PHASE -TRANSFER  PROCESSES  AND  PROPERTIES 

GROUP  11.  INITIALIZATION  OF  VARIABLE  OR  POROSITY  FIELDS 
FIINIT( HI  )=HAMB 

FIINIT(TMP1)=TAMB  J  FIINIT( DENl )=1 .0 
FIINIT(KE)=KEINA  }  FIINITl EP J=EPINA 
FIINIT1P1)=0.0 
FIINIT(W1)=0.0 
FIINIT(VU=0.0 
FIINIT(U1)=0.0 

PATCH1WINIT,INIVAL,1,5,1,8,1S,30,1,1) 

INIT(WINIT,W1,0.0,0.2»WJET) 

INIT( WINIT ,H1 ,0 . 0 ,CSUBP*0 . 2«TJET ) 

INIT( HINIT ,KE ,0 . 0 ,KEINJ*0 . 2 ) 

INIT(WINIT,EP,0.0,EPINJ»0.2) 
»«*  SEE  GROUP  13  FOR  POROSITY  FIELDS 

GROUP  12.  CONVECTION  AND  DIFFUSION  ADJUSTMENTS 

GROUP  13.  BOUNDARY  CONDITIONS  AND  SPECIAL  SOURCES 
«*DEFINE  SOURCES  AND  SINKS  FOR  JET  ENGINE** 
CONPOR(0.0,CELL,1,1,<+,S,12,1^) 
PATCH(  JETIN, HIGH,  1,1, '♦.5,11, 11  »1,1) 
COVAL( JETIN,P1,FIXFLU,-MJET) 
COVAL( JETIN,H1,0NLYMS,SAME ) 
PATCHl  JET0UT,L0W,1,1,<»,5,15,15,1,1) 
COVAL( JET0UT,P1,FIXFLU,MJET) 
COVALI JET0UT,H1,0NLYMS,HJET) 
C0VAL( JETCUT ,W1 ,ONLYMS,WJET ) 
COVAL( JETOUT ,KE ,ONLYMS ,KE INJ ) 
COVAL(JETOUT,EP,ONLYMS,EPINJ) 

«*EAST  BOUNDARY  TO  ATMOSPHERE** 
PATCH!  OUTSIDE, EAST,  10, 10,1, 12, 1,<»0, 1,1) 
COVAL( OUTSIDE, PI, 0.1, 0.0) 
COVAL( OUTSIDE , HI , FIXVAL ,HAM8 ) 

*«TOP  OF  HUSH  HOUSE  AND  TUBE  INCLUDING  TUBE  LIP** 
CONPORi 0.0, NORTH, 1,9, 12, 12, 1,19) 
C0MP0R(0.0,N0RTH,l,5,-8,-8,18,27) 
CONPORI 0.0, NORTH, 1,5, -2, -2, 35, 35) 
C0NPOR(  0.0, NORTH,  1,5, -<♦,-<♦,  36, 36) 
CONPORI 0.0, NORTH, 1,5, -5, -5, 37, 37) 
CONPOR(0.0,NORTH,l,5,-6,-6,38,38) 
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CONPORl0.0,NORTH,l,S,-7,-7,39,59) 

♦HtTOP  OF  MODEL  TO  ATMOSPHERE** 
PATCH(TOPOaT,NORTH,l,10,12.12.20,*0,l,l) 
COVAL(TOPOUT,P1,0.1,0.0) 
COVAL( T0P0UT,H1 ,ONLYMS,HAMB ) 
PATCH(TOPP,NORTH,10,10,12,12,1»19,1,1) 
COVAL(TOPP,P1,0.1,0.0) 
C0VAL(T0PP,H1,0NLYMS,HAMB) 

w*FLOOR  OF  HUSH  HOUSE.  TIAE  AND  GROUND** 
CONPOR(0.0»SOUTH,l,10,l,-l,I,<iO) 

»*FORE  AND  AFT  BULKHEADS  IN  HUSH  HOUSE  AND  TUBE* 
CONPOR( 0.0, LOW, 1,9, 1,12, -1,-1) 
CCNPOR(0.0,LOH,6,9,1,12,-20,-20) 
CONPORI 0.0, LOW, 1,5,9,12,-20, -20) 
CONPOR(0.0,LOW,l,5,l,2,-3S,-5S) 
CONPOR(0.0,LOW,l,5,5,'f,-36,-36) 
CONPOR(0.0,LOH,l,5,S,5,-37,-37) 
CONPOR(0.0,LOW,l,5,6,6,-S8,-58) 
CONPOR(0.0,LOW,l,5,7,7,-39,-39) 
CONPOR(0.0,HIGH,l,S,8,8,-39,-39) 

**VERY  FRONT  MODEL  STRIP  TO  ATMOSPHERE** 
PATCH( FRONT, LOW, 10, 10, 1,12, 1,1, 1,1) 
COVALl FRONT, PI, 0.1, 0.0) 
COVAH FRONT, HI, ONLYMS,HAMB) 

«*THE  VERY  END  OF  OUR  MODEL  TO  ATMOSPHERE** 
PATCH(ENDEN0, HIGH,  1,10,1, 12, <»0,<»0, 1,1) 
COVAL(ENDEND,P1,0.1,0.0) 
C0VAL(ENDEND,H1,0NLYMS,HAMB) 

**OUTSIDE  EAST  HALL  IN  HUSH  HOUSE** 
CONPOR(0.0,EAST,-9,-9,9,12,l,19) 
CONPOR(0.0,EAST,-9,-9,l,8,l,l) 
CONPORl0.0,EAST,-9,-9,l,8,3,3) 
CONPORI0.0,EAST,-9,-9,l,8,5,8) 
CONPOR(0.0,EA3T,-9,-9,l,8,10,ll) 
CONPOR(0.0,EAST,-9,-9,l,8,13,14) 
CONPOR(0.0,EAST,-9,-9,l,8,18,19) 

*«OUTSIDE  OF  TUBE  PLUS  TUBE  LIP  IN  HUSH  HOUSE** 
CONPOR(0.0,EAST,-5,-5,l,8,18,39) 

GROUP  l<t.  DOWNSTREAM  PRESSURE  FOR  PARAB=.TRUE. 

GROUP  IS.  TERMINATION  OF  SWEEPS 
FSWE£P=li  LSWEEP=180 

GROUP  16.  TERMINATION  OF  ITERATIONS 

GROUP  17.  UNOER-RELAXATION  DEVICES 
REAHD£LT,FACT) 
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0ELT=50 . 0/i  NZ»0 . 2»«WJET ) 

FACT=10. 

FACT=2.5 
FACT=0.7S 

FACT=3.S 

FACT=S. 

FACT=0.50 

FACT=0.2S 
RELAX(P1,LINRLX,0.3) 
RELAX! Ul ,FALSDT,FACT»DELT ) 
RELAXIV1,FALSDT,FACT«DELT) 
RELAX! Wl ,FALSDT,FACT«OELT ) 
RELAX! H1,FALSDT,2.S«FACT«DELT) 
RELAX! KE,FALSDT,FACT*DELT) 
RELAX! EP,FALSDT,FACT«OELT) 

GROUP  18.  LIMITS  ON  VARIABLES  OR  INCREMENTS  TO  THEM 
VARMINIUl)=-200. 
VAW1AXIU1)  =  1000. 
VARMIN(Vl)=-300. 
VARMAXIV1)=1000. 
VARMINIW1)=-100. 
VARMAX1W1)=1S00. 
VARMIN!TMP1)=2S0 
VARMAX!TMP1)=TJET 
VARMIN(H1)=1000. 
VARMAXI HI  )=CSUBP»TJET 
VARMINIPl)=-3.5E+05 
VARMAXI PI  )=  3.5E+0S 
VARMIN(RH01)=.OS 
VARMAXI RHOl) =5. 

VARMINIKE)=l.E-08 

VARMAX(KE )=1.3*KEJET 

VARMINIEP)=l.E-08 

VARMINIVIST)=0.1 

VARMAXIEP)=1.3»1KEJET»«1.S) 

GROUP  19.  DATA  COMMUNICATED  BY  SATELLITE  TO  GROUND 

GROUP  20.  PRELIMINARY  PRINT -OUT 

GROUP  21.  PRINT -OUT  OF  VARIABLES 

GROUP  22.  SPOT-VALUE  PRINT -OUT 

»«MONITOR  THE  OUTPUT  AT  DESIGNATED  POINTS** 
IXM0N=1  i    IYM0N=4  i  I2M0N=16 
TSTSNP=1  i    LUPR3=6 
ABSI2=1.0  i    0RSIZ=1.0 

GROUP  23.  FIELD  PRINT-OUT  AND  PLOT  CONTROL 
NUMCLS=10*HXPRIN 
ITAEL=3  i    NPLT=1  )  N2PRIN=1 
IPLTL=LSHEEP 
NPRINT=3 
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GROUP  Z^.  DUMPS  FOR  RESTARTS 
SAVE=T 
RESTRT(ALL) 
STOP 
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APPENDIX    B  

SEH-JECTED    PHOEITICS    OUTPUT 


WMmmMMM-MROItHMWMMHiOlllllilMMWWWW* 


CCCC   HHH 

CCCCCCCC    HHHHH 

CCCCCCC    HHHHHHHHHH 

CCCCCCC   HHHHHHHHHHHH 

CCCCCC   HHHHHHHHHHHHH 

CCCCCCC   HHHHHHHHHHHH 

CCCCCCC   HHHHHHHHHH 

CCCCCCCC   HHHHH 

CCCC   HHH 


PHOENICS     VERSION   1.3>    03  SEPT   1986 

(C)    COPYRIGHT    leS^t 

CONCENTRATION   HEAT   AND   MOMENTUM   LTD 

ALL   RIGHTS   RESERVED. 

CHAM   LTD,   BAKERY   HOUSE,   40   HIGH  ST 

WIMBLEDON,    LONDON,   SW19  SAU 

TEL:    01-9<+7-76Sl>   TELEX:    928S17 

FACSIMILE:    01-879-2^97 

THE   OPTION   LEVEL   IS        -18 


THIS  CODE   MAY   ONLY   BE   USED   UNDER   THE   TERMS  AND  CONDITIONS 
OF  A   LICENCE   AGREEMENT  WITH   CHAM   LTD. 

REPLICATION  OF   THIS  CODE    IS  PROHIBITED  UNLESS 


SPECIFICALLY   AUTHORISED   IN  WRITING  BY   CHAM   LTD. 

GREXl  OF    15/07/86   HAS   BEEN   CALLED 

FORMATTED   SATLIT   DATA   READ    FROM   DFIO    FOR   IRUN=      1 

*«* STORAGE    INFORMATION «**  ' 

F   DIMNSN=200000   OCCUPIED=12S707   ESTIMATED   MINIMUM   DIMNSN=   23267 
DEP  VRBL=  67200   OLD  VRSL=  0   3D   COEFF=   23880   SD   DVDPS=   l<t<tOO 

INITIAL   FIELDS   READ   FROM   DF09 

GROUP   1.    RUN  TITLE    S   NUtlBER 


TEXT(NAS  JACKSONVILLE   HUSH   CHAMBER 


) 


fr  »Hr ft  H  A  N  H  ^"^  ^w"W^t^r^c  W^UHmTt  >i 


STEADY 


CARTES 

n;< 

XULAST 


IRUNN        =  1 

Vtl(#Wl(MlVRtt9lft)(viA9t  tMt  ?HH<  ft  Jt  VV  V 

GROUP   2.    TRANSIENCE 

T 

W  It  It  K  n  ^  ^  W^W  w  W  W  B  ^  ^  W  M  n  H  M  R  n  H  W^V  n  n  H  K  K  H  ■rWT^Tnr  K  H  wt^f^rtm 

GROUP   3.    X-DIRECTION   GRID   SPACING 
T 
10 
l.OOOE+00 

METHOD   OF    PAIRS   USED    FOR   GRID  SETTING 

XFRAC    (       1)    =    -l.OOOE+00    jXFRAC  I       2)    = 

XFRAC    (       3)    =      <».000E*00    jXFRAC  (       <♦ )    = 

XFRAC    (      5)    =      3.000E+00    jXFRAC  (       6)    = 

XFRAC    (       7)    =      2.000E*00    >XFRAC  (       8)    = 

GROUP  4.    Y-DIRECTION  GRID   SPACING 
NY  =  12 

YVLAST      =   1. CODE ♦GO 

METHOD   OF   PAIRS   UGED    FOR   GRID   SETTING 


3.100E-01 
5.<t7SE-01 
2.000E+00 

mooE+oo 
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YFRAC  (  1)  = 

YFRAC  (  3)  = 

YFRAC  (  5)  = 

YFRAC  (  7)  = 


-3.000E+00  )YFRAC  ( 

2.000E+00  jYFRAC  ( 

3.000E+00  JYFRAC  ( 

<».OOOEtOO  »YFRAC  J 


2) 

s 

5 

620E-01 

<♦) 

= 

3 

l<tOE-01 

6) 

s 

5 

620E-01 

8) 

= 

7 

500E-01 

^^^  n  M  K  MRMMMHRMMnX  W  ft  ^r^^  W  M  H  H  M  V 

GROUP  S.    Z-DIRECTION  GRID   SPACING 
PARAB        =  F 

NZ  =40 

2WLAST      =    l.OOOE+00 

METHOD   OF   PAIRS  USED   FOR  GRID  SETTING 


2FRAC    ( 

1) 

= 

-1 

000E*00 

>ZFRAC    ( 

2) 

r 

1 

OOOE+00 

2FRAC 

3) 

s 

1 

OOOE+00 

jZFRAC    ( 

<+) 

s 

1 

700E+00 

2FRAC    ( 

S) 

= 

1 

OOOE+00 

)ZFRAC    ( 

6) 

= 

2 

OOOE+00 

2FRAC 

7J 

= 

1 

OOOE+00 

jZFRAC    ( 

8) 

r 

1 

700E+00 

ZFRAC 

9) 

= 

2 

OOOE+00 

jZFRAC    ( 

10) 

= 

3 

OOOE+00 

ZFRAC 

11) 

= 

1 

OOOEtOO 

JZFRAC    ( 

12) 

= 

1 

500E+00 

ZFRAC    ( 

13) 

= 

1 

OOOE+00 

jZFRAC    ( 

1<») 

= 

1 

OOOE+00 

ZFRAC 

IS) 

= 

1 

OOOE+00 

jZFRAC    ( 

16) 

= 

1 

700E+00 

ZFRAC 

17) 

= 

2 

OOOE+00 

jZFRAC    J 

18) 

= 

1 

OOOE+CO 

ZFRAC    ( 

19) 

= 

1 

OOOE+00 

iZFRAC    ( 

20) 

= 

1 

700E+00 

ZFRAC    ( 

21) 

= 

2 

OOOE+00 

jZFRAC    ( 

22) 

: 

1 

OOOE+OO 

ZFRAC    ( 

23) 

= 

3 

OOOE+00 

jZFRAC    ( 

24) 

= 

S 

667E-01 

ZFRAC    ( 

25) 

= 

2 

OOOE+00 

JZFRAC    ( 

26) 

= 

3 

500E-01 

ZFRAC    ( 

27) 

= 

7 

OOOE+00 

jZFRAC    1 

28) 

= 

3 

OOOE+00 

ZFRAC    ( 

29) 

= 

1 

OOOE+00 

jZFRAC    ( 

30) 

- 

2 

OOOE+OO 

ZFRAC    1 

31) 

s 

<♦ 

OOOE+00 

iZFRAC    ( 

32) 

- 

7 

SOOE-01 

ZFRAC    ( 

33) 

= 

<► 

OOOE+00 

jZFRAC    ( 

34) 

= 

5 

7S0E-01 

ZFRAC    ( 

35) 

= 

<♦ 

OOOE+00 

jZFRAC    ( 

36) 

= 

£ 

OOOE-01 

ZFRAC    I 

37) 

= 

1 

OOOE+00 

)ZFRAC    { 

38) 

s 

3 

OOOE+00 

w  w  ^^  ^  w  w  M  Ann  w^  ^r^r^T^r^^^  ^  ^^  m.  w  ^^^  nnnnwnMMV 

GROUP    6.    BFC    AND   GEOMETRY   CHA^JGES 

GROUP    7.    VARIABLES:    STORED .SOLVED , NAMED 
ONEPHS      =  T 

NAME(  1)  =P1 
NAME(  5)  =V1 
NAME  I  12)  sKE 
NA>1E(14)  =H1 
NAME  120)  : 
NAME (22)    : 


Ht  w  w  H  jl  W"W^T  ft  ft  tt  w 


jNAME(     3)  =U1 

>NAME(     7)  =W1 

>NAME( 13  >  =EP 

>NAM£( 19)  =HPOR 

NPCR    ;rJAME(  21)  =EPOR 

VPOR    »NAME123)  =VIST 


NAME(2<+)  =RH01  )NAME!25)  =TMP1 

SOLUTNIPl   ,Y,Y,Y,N,N,N) 

SOLUTNCUl   ,Y,Y,N,N,N,N) 

SOLUTNIVl   ,Y,Y,N,N,N,N) 

SOLUTNIHI   ,Y,Y,N,N,N,N) 

SOLUTN(KE   ,Y,Y,N,N,N,N) 

SOLUTN(EP   ,Y,Y,N,N,N,N) 

SOLUTN(HI   ,Y,Y,N,N,N,N) 

SOLUTN(KPOR,Y,N,N,N,N,N) 

SOLUTN(NPOR,Y,N,N,N,N,N) 

SOLUTN(EPOR,Y,N,N,N,N,N) 

SOLUTN(VPOR,Y,N,N,N,N,N) 

SOLUTN(VIST,Y,N,N,N,N,N) 

S0LUTN(RH01,Y,N,N,N,N,N) 

S0LUTN(TMP1,Y,N,N,N,N,N) 

DENl     =     Z(* 

VIST    =     23 

EPOR     =     21  )HPOR    =     19  jNPOR 

TEMPI    =     25 

K^nnAHH^  H  n  H  H  M  A  ^^€  R  n  M  ^^  Tf 

GROUP  8.  TERMS   S  DEVICES 


20  JVPOR 


22 


TERMS  (PI 

,Y,Y,Y,N,Y,N) 

TERMS  (Ul 

,Y,Y,Y,N,Y,N) 

TERMS  (VI 

,Y,Y,Y,N,Y,N) 

TERMS  (HI 

,Y,Y,Y,N,Y,N) 

TERMS  (KE 

,N,Y,Y,N,Y,N) 

TERMS  (EP 

,N,Y,Y,N,Y,N) 

TERMS  (HI 

,Y,Y,Y,N,Y,N) 

DIFCUT   = 

S.OOOE-01  iZDIFAC 

=  1 

OOOEtOO 

GALA     = 

F   lADDDIF  = 

F 

NEWRHl   = 

T 

NEHENT   = 

T 

«»-)t»*»»»<t»»»«»<HHt»»<t  K  KM  M  »  mH 

****^ 

(«»«-»«  «-K- 

GROUP  9. 

PROPERTIES 

RHOl     =- 

-1.016E+0'*  jTMPl 

--1 

013E+0'*  >EL1 

=  -1 

015E*0'+ 

RHOIA    = 

O.OOOE+00  JRHOIB 

-  Z 

<>96E-03  >RH01C 

=  7 

1<»3E-01 

PRESSO   = 

1.013E+C5 

TMPIA    = 

l.OOOE-20  )TM?1B 

=  9 

960E-0<»  jTMPlC 

=  0 

OOOE'^00 

TEMPO    = 

O.OOOE+00 

ELIA     = 

O.OOOE+OO  ,EL1B 

=  0 

OOOE+OO  )EL1C 

=  0 

OOOE+OO 
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ENUL  =   l.OOOE-05    lENUT 

ENUTA        =   O.OOOE+OO    jENUTB 
DRHIDP      =-1.016E+0^ 


=  -1.01'»E*0«* 

s   O.OOOEtOO    jENUTC 


O.OOOE^OO 


l.OOOE+00  iPRNDTLtVl 
l.OOOE+00  iPRNOTLIKE 
l.OOOE+OO  )PRNDTL(H1 
l.OOOE+00  jPRT  IVl 
l.OOCE+00  >PRT  (KE 
l.Sl'tE+OO    )PRT         (HI 

KM  MMMM^KHKH  ^TT*  W^WTnr  sT  W^  W  W  R  W  ^  W  W^W^  M  R  M  M  ^  M  M  R  M  H  W  M  ^t^K  ^  H  H  H  A  H  ^TtT^V^  M  R  ^FW^r^^WyK 

GROUP   10. INTER-PHASE   TRANSFER   PROCESSES 

GROUP   11. INITIALIZE   VAR/POROSITY    FIELDS 

FIINITIPl      )    =    -1.022E+0<*  jFIINITIUl       ) 

-1.022E*0^  jFIINIT(Wl       ) 

-1.022E+0'*  >FIINIT(EP       ) 

-1.022E*0<f  iFIir<IT(HPOR) 

-1.022E+04  jFIINIT(EPOR) 

-1.022E*0<*  JFIINITIVIST) 

-1.022E+0'+  iFIINIT(TMPlJ 


PRNDTLIUl 
PRTJOTLIWI 
PRNDTLIEP 
PRT  (Ul 
PRT  (Ml 
PRT        (EP 


FIINITIVI  ) 
FIINITIKE  J 
FIINIT(H1  ) 
FIINIT(NPOR) 
FIINIT(VPOR) 
FIINIT(RHOl) 


l,000E+00 
l.OOOEtOO 
l.OOOE+00 
l.OOOE+00 
1,000E*00 
l.OOOE+00 


-1.022E+04 
-1.022E+0<+ 
-1.022E+0'+ 
-1.022E+0^ 
-1.022E+0'+ 
-1.022E+0'+ 
-1.022E+0<+ 


PATCH(S0UTH1       .INIVAL,         1,         1,         3,         3,      12,      1^ , 
INITISOUTHI      ,V1      ,    O.OOOE+00,    O.OOOE+00) 
INIT(S0UTH1      ,NPOR,    O.OOOE+00,    O.OOOE+00) 


1>        1) 


PATCH(L0W1  ,INIVAL,         1,  1,         4,         5, 

INIT(L0N1  ,H1       ,    O.OOOE+00,    O.OOOE+OO) 

INIT(L0W1  ,HPOR,    O.OOOE+OO,    O.OOOE+OO) 


11,      11,        1, 


1) 


PATCH 
INITl 
INITI 
INITl 
INITl 
INITl 
INITl 
INITl 
INITl 
INITl 
INITl 


ICEL 
CELL 
CELL 
CELL 
CELL 
CELL 
CELL 
CELL 
CELL 
CELL 
CELL 


,INIVAL,         1,         1,         <f,         5,      12>      1<+,         1,         1) 

,P1      ,    O.OOOE+00,  O.OOOE+OO) 

,U1      ,    O.OOOE+00,  O.OOOE+OO) 

,V1      ,    O.OOOE+OO,  O.OOOE+OO) 

,W1       ,    O.OOOE+OO,  O.OOOE+OO) 

,KE      ,    O.OOOE+OO,  O.OOOE+OO) 

,EP      ,    O.OOOE+00,  O.OOOE+OO) 

,H1       ,    O.OOOE+00,  O.OOOE+00) 

,HP0R,    O.OOOE+OO,  O.OOOE+OO) 

,NP0R,    O.OOOE+00,  O.OOOE+00) 

,EPOR,    O.OOOE+00,  O.OOOE+CO) 
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PATCH(EAST21      ,INIVAL,         9,        9,        1,        8, 
INITIEAST21      ,U1      ,    O.OOOE+00,    O.OOOE+00) 
INITIEAST21      >EPOR,    O.OOOEi'OO,    O.OOOE'^00) 


3,        i,        1,        I) 


PATCH! EAST22      .INIVAL,         9,        9,        1,        8, 
INIT(EAST22      ,U1      ,    O.OOOE+00,    O.OOOE+OOJ 
INIT(EAST22      ,EPOR,    O.OOOE^OO,    O.OOOEtOO) 


S, 


1,        1) 


PATCH(EAST23      .INIVAL,         9,        9,        1,        8, 
IMIT(EAST23      ,U1      ,    O.OOOE+00,    O.OOOE+00) 
INITIEAST23      ,EPOR,    O.OOOE+00,    O.OOOE+00) 


10,      11,        1, 


1) 


PATCH(EAST2<>       , INIVAL,         9,         9,         1,         8,      13,      1<», 
INIT(EAST2<»      ,U1       ,    O.OOOE+OO,    O.OOOE+OO) 
INIT(EAST24      ,EPOR,    O.OOOE+OO,    O.OOOE+00) 


1, 


IJ 


PATCH(EAST2S      , INIVAL,         9,        9,        1,        8,      18, 
INIT(EAST25      ,U1       ,    O.OOOE+00,    O.OOOE+00) 
INIT(EAST25      ,EPOR,    O.OOOE+00,    O.OOOE+00) 


19,        1, 


1) 


PATCH! EAST26      .INIVAL,        5,        S,        1,        8,      18,      39, 
INIT(EAST26      ,U1      ,    O.OOOE+OO,    O.OOOE+00) 
INIT(EAST26      ,EPOR,    O.OOOE+OO,    O.OOOE+OO) 
INIADD      =  F 

LUFI  =  9 

NAMFI        =DF09 

iT^n  K  H  n  yr^f  M  w  ^  ^r^T  n  n  M  ■  wvf  Tucw^TnT^r^  n  m  h  n  n  ^rw^r^w^ndi 

GROUP   13.    BOUNDARY    i   SPECIAL  SOURCES 


1,        1) 


PATCHl RESOURCE, PHASEM,  1,  10,  1,  12, 
COVAL ( KESOURCE ,KE  , -1 . OlSE  +0^ ,  -1 . OlSE  +0^ ) 
COVALlKESOURCE,EP      ,-1.015E+0't,-1.01SE+0'+) 


1,      ^0,        1, 


1) 


PATCH!  JETIN         ,HIGH      ,         1,        1,        <♦,        S,      11, 
COVALIJETIN        ,P1      ,    1.000E-10,-1.977E+02) 
COVALIJETIN         ,H1      ,    O.O0OE+O0,-1.026E+O4 ) 


11,        1, 


1) 


PATCH! JETOUT       ,LOW         ,         1,         1,        ^,         S,      15 j 

COVAL! JETOUT      ,P1      ,    I.OOOE-IO,    1.977E+02) 


IS, 


1) 
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C0VAL(HWALH2  ,U1 
C0VAL(HWALL12  ,V1 
CCVAL(HWALL12  ,KE 
C0VAL(HWALL12    ,EP 


,-1.01lB*0it,  O.OOOE^OO) 
,-1.013E+0'+,  O.OOOEtOO) 
,-1.01SE+0^,-1.013E+0<»J 
j-l.OlSE+O-^j-l.OlJE+O^) 


20>      20,        1,        1) 


3<»,      3^,        1, 


35,      35,        1, 


1) 


1) 


PATCH! LKALL12  ,LWALL    ,         1,         5,         9,       12, 

C0VAL(LWALL12  ,U1       ,-l.OlSE+O^,    O.OOOE+00) 

C0VALILWALL12  ,V1      ,-1.013E+0^,    O.OOOE+00) 

C0VAL(LWALL12  ,KE       ,-1 .013E+0^ ,-1 , 013E*0^  ) 

C0VAL(LWALL12  ,EP      ,-1.013E*04,-1.013E*0<» ) 

PATCH(HWALL13  ,HWALL    ,         1,         S,         1,         2, 

C0VAL(HWALL13  ,U1      ,-1.013E*0<*,    O.OOOE+00) 

C0VAL(HWALL13  ,V1      ,-1 .013E tO<+,    O.OOOE^OO) 

C0VALIHWALL13  ,KE       ,-1 .013E  +  0<+,-l .  013E*0<+ ) 

C0VAL(HHALL13  ,EP      ,-1 .013E*04,-1.013E*0^  ) 

PATCH(LKALL13  ,LWALL    ,        1,        5,        1,         2, 

C0VAL(LWALL13  ,U1      ,-1 .013E+0<i,    O.OOOE+00) 

C0VAL(LWALL13  ,V1       ,-1 .  OlSE+O**  ,    O.OOOE+00) 

C0VAL(LWALL13  ,KE       ,-1. 013E*0^,-1.013E*0<*  ) 

COVAKLWALLIS  ,EP      ,-1.013E tO<»,-l .013E+0'* ) 

PATCHIHWALLI^  ,HWALL    ,         1,         5,         3,         <♦ , 

COVALIHWALU^  ,U1       ,-1 .013E+0'^  ,    O.OOOE+OO) 

COVAKHl'tALLl'*  ,V1       ,-1 .013E*0'> ,    O.OOOE+00) 

COVA!.(HKALLl't  ,KE       ,-1 .013E+0'» , -1 .015E+0t ) 

C0VAL(H)-<ALL1<+  ,EP       ,-1 .013E+0<^,-1.013E+0<»  ) 


PATCH!  LWALLl^  ,LHALL    ,         1,         5,         3,         <* ,      36,      36,         1,         1) 

COVALILKALLI^  ,U1      ,-1.013E+0^,    O.OOOE+OOI 

COVALILWALLI'*  ,V1      ,-l,013E+0^,    O.OOOE+OO) 

C0VAL(LWALL1<+  ,KE       ,-1 .013E+0^  ,-1 .015E  +  0<+ ) 

C0VAL(LKALL1'+  ,EP       ,-1 .013E+0^,-l  .013E+0^  ) 

PATCHIHWALLIS  ,HWALL    ,        1,        5,        S,        S,      36,      36,        1,        1) 

COVALtHWALLlS  ,U1      ,-1.013E+0<* ,    O.OOOE+OO) 

COVALIHVJALLIS  ,V1      ,-1.013E+0^,    O.OOOE+OO) 

COVALCHi^ALLlS  ,KE       ,-1 .  OlSc+O'*,-!  .013E +0-+ ) 

COVAKW'^ALLIS  ,EP      ,-1. OlSE+O't ,-l.C13E+0'+ ) 


35,      35,        1,        1) 
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PATCHfLWALLlS 
C0VALILWALL15 
COVALILHALLIS 
COVAK  LWALLIS 
COVALC  LHALL15 


,LWALL 

,U1 

.VI 

»KE 

,EP 


PATCH{HWALL16  ,HWALL 

COVALC HWALL16  ,U1 

C0VALIHWALL16  ,V1 

C0VAL1HWALL16  ,KE 

C0VAL(HHALL16  ,EP 

PATCH(LWALL16  ,LWALL 

C0VAL(LWALL16  ,U1 

C0VAL(LHALL16  ,V1 

C0VAL(LWALL16  ,KE 

COVALC LWALL16  ,EP 

PATCH(HWALL17  ,HWALL 

C0VAL(HWALL17  ,U1 

C0VAL(HWALL17  ,V1 

C0VAL(HWALL17  ,KE 

C0VALIHWALL17  ,EP 

PATCH! LWALL 17  ,LWALL 

C0VAHLWALL17  ,U1 

C0VALILWALL17  ,V1 

C0VAL(LWALL17  ,KE 

C0VAL(LHALL17  ,EP 

PATCH(HWALL18  ,HWALL 

C0VAL(HWALL13  ,U1 

C0VAL(HWALL18  ,V1 

C0VAL(HKALL18  ,KE 

COVAH  HWALL18  ,EP 

PATCHI LWALLIS  ,LHALL 

C0VAL( LWALLIS  ,U1 

COVAH  LVULLIS  ,V1 


1,  5>  S,  S, 
-1.013EtO<f,  O.OOOE^OO) 
-1.013E«04,  O.OOOEfOO) 
-l.OlSE+C+.-l.OlSE+O^t) 
-1.013E*0';,-1.012E*0<») 

>  1»        S>        o»         6> 

-1.015E*0^,  O.OOOE+00) 
-1.013E+0<+,  O.OOOE+00) 
-1.013E+0<»,-1.013E+0^) 
-1.013E*0<*,-1.013E*0't) 

>  1>    5>    oj    o» 

-1.013E+0'*,  O.OOOE+00) 
-LOISE+C*,  O.OOOE+00) 
-1.013E*04,-1.015E+0'+) 
-1.013E*0'»,-1.013E*0'+) 

,  1,  S,  7,  7, 
-1.013E*0<f,  O.OOOEfOO) 
-1.013E+0«+,  O.OOOE+OO) 
-1.013E+0^,-1.013E+0^) 
-1 .  013E  tCt , -1 .  013E+0<»  ) 

1,  S,  7,  7, 
-1.012E*0^,  O.OOOE+00) 
-1.013E+0^,  O.OOOE+00) 
-1.013E*0^,-1.013E+0^) 
-1.013E*04,-1.015E*04) 

»   1>   S>   8»    8» 

-1.013E+0'*,  O.OOOE+00) 
-1.013E+0^,  O.OOOE+00) 
-1.013Et0^,-1.013E+0'f) 
-1.013E*0^,-1.013E*0'*) 

»        1»        S»        8>        8> 

-1.013E*0^,  0.000E*00> 
-1.013E+0''>,    O.OOOE+00) 


37,      37,        1, 


1) 


37,      37, 


1,        1) 


38,      38,        1,        1) 


38,      38,        1, 


39,      39,        1, 


39,      39,        1, 


<+0,      <»0,        1, 


1) 


1) 


1) 


1) 
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COVALI LWALL18    ,KE 
C0VAL(LWALL18    ,EP 


,-1.013E*0<»,-1.013E+0<>) 
,-1.013E*0<»,-1.013E*0't) 


PATCH( FRONT  ,LOW         ,      10,      10,         1,      12, 

COVAL( FRONT  ,P1      ,    l.OOOE-01,    O.OOOE+00) 

COVALI  FRONT  ,H1      ,    O.OOOE^'OO,    2.962E^05) 

PATCH! ENDENO  ,HIGH      ,        1,      10,        1,      12, 

COVAKENDENO  ,P1      ,    l.OOOE-01,    O.OOOE+00) 

COVAKENDENO  ,H1      ,    O.OOOE  +  00,    2.962E+05) 

PATCH(EWALL19  ,EWALL    ,         9,         9,         9,      12, 

C0VAL(EWALL19  ,V1       ,-1.013E*04,    O.OOOE+00) 

C0VAL(EHALH9  ,W1       ,-1.013E+0^,    O.OOOE  +  00) 

C0VAL(EWALL19  ,KE       ,-1 .  OlSE+O'*, -1 .  013E+0'*  ) 

C0VAL(EWALL19  ,EP      ,-1 .013E  fO^, -1 .013E+0<*  ) 

PATCH(KHALL19  ,J<WALL    ,      10,      10,         9,      12, 

C0VAL(KWALL19  ,V1       ,-1 .  013E*0<*  ,    O.OOOE+00) 

C0VAL1WHALL19  ,W1       ,-1 .013E*0<+,    O.OOOE+OC) 

C0VAL(HWALH9  ,KE       ,-1 .  013E*0h,-1.013E*0'*  ) 

C0VAL1WJ-1ALL19  ,EP       ,-1 .  013E*0^,-1 .013E  +  0<»  ) 

PATCH(EWALL20  ,EWALL    ,        9,        9,        1,        8, 

COVAKEWALLZO  ,V1      ,-1.013E+04,    O.OOOE+00) 

COVAL(EWALL20  ,W1       ,-1 . OlSE+CV,    O.OOOE+00) 

COVAL(EWALL20  ,KE       ,-1 .013EfO<+,-l .  015E+0<+ ) 

COVALIENALL20  ,EP       ,-1 .  OlSE+O"^,-!  .013E  +  0^  ) 

PATCH(WWALL20  ,WWALL    ,      10,      10,        1,        8, 

COVAL(WWALL20  ,V1       ,-1.013E+04,    O.OOOE+00) 

C0VAL(WWALL20  ,W1      ,-1  .OliE+O'f ,    O.OOOE+00) 

COVAL(WHALL20  ,KE       ,-1 .  013E+0'*  ,-1 .  013E  +  0^  ) 

COVAL(WHALL20  ,EP       ,-1 .  013E+0't, -1 .  013E*0'+ ) 

PATCH(EHALL21  ,ENALL    ,         9,         9,         1,         8, 

C0VAL(EKALL21  ,V1       ,-1.013E  +  0'» ,    O.OOOE  +  00) 

C0VAL(EWALL21  ,W1       ,-1  .OlSE+O"* ,    O.OOOE  +  00) 

COVAHEHALLZl  ,KE       ,-1 .  013E+0<»  ,-1 .013E*0<>  I 

C0VALIENALL21  ,EP       ,-1 .  013£+0<^,-l .  013E*0^  ) 


1>         1, 


<tO,     w, 


1) 


1,     1) 


1,      19,        1, 


1,      19,        1, 


1,        1,        1, 


1.        1,        1, 


3,        3,        1, 


1) 


1) 


1) 


1) 


1) 
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PATCH(WWALL21  ,WWALL    ,      10,      10,         1,        8, 

C0VAL(WWALL21  ,V1      ,-1 . 013E*0<»,    O.OOOE+00) 

C0VAHWWALL21  ,H1      ,-1. 013E*0<f ,    O.OOOE^OO) 

C0VAL(WWALL21  ,KE       ,-1 .OlSE+C*,-! .013E+04 J 

C0VAL(WWALL21  ,EP      ,-1.  013E*0<*,-1.013E*0<i  ) 


3,   3,   1, 


1) 


PATCH! EWALL22  ,EWALL  ,  9,  9,  1,  8, 
COVALIEWALL22  ,V1  ,-1.013E+0^,  O.OOOE+00) 
COVAL(EWALL22  ,W1  ,-1.013E*0^,  O.O00E*O0) 
COVAL(EHALL22  ,KE  ,-1 .013E  +  0^,-1 .013E+0'+ ) 
C0VAL(EWALL22  ,EP   ,-1.013E*0^,-1.013E*0'* ) 


1) 


PATCH(WWALL22  ,WWALL  ,   10,   10,    1,    8, 

COVAL(WWALL22  ,V1   ,-1.013E*04,  O.OOOE+00) 

COVAL(HWALL22  ,W1   ,-1 .  OlSE+O"^,  O.OOOE+00) 

COVAL(WWALL22  ,KE   ,-1 .  013E*0<»,-1 .013E*0<+ ) 

C0VALIKWALL22  ,EP   ,-1.013E*0^,-1.013E+0'* ) 


5,    8,    1,    1) 


PATCH(EWALL23  ,EWALL  ,    9,    9,    1,    8, 

COVAL(EWALL25  ,V1   ,-1. 013E+0<*,  O.OOOEtOO) 

C0VAL(EWALL23  ,W1   ,-1 . 013EtO<i,  O.OOOE  +  OOJ 

COVALIEWALL23  ,KE   ,-I.013E+0't,-l  .013E  +  0^  ) 

COVALIEWALL23  ,EP   ,-1.013E+0'*,-1.013E*0<t ) 


10,   11,    1, 


1) 


PATCH(WWALL23  ,WWALL  ,   10,   10,    1,    8, 

C0VAL(WWALL23  ,V1   ,-1.013E+0^,  O.OOOEtOO) 

COVALIWWALL23  ,W1   ,-1.013E+0^,  O.OOOE+OO) 

C0VAL(WWALL23  ,KE   ,-1 .  013E  +  04  ,-1 .013E*0'*  ) 

C0VAL(W;'1ALL23  ,EP   ,-1 .013E*0^,-1 .013E*0^  ) 


10,      11, 


1,        1) 


PATCH!  EWALL2<»  ,EHALL    ,         9,        9,         1,        8, 

C0VALIEWALL2^  ,V1       ,-1.013E+0^,    O.OOOE+00) 

C0VALIEWALL2<+  ,H1       ,-1.013EtO^,    0.000E*00) 

C0VAL!EWALL2'+  ,KE       ,-1 .  OlSE+O'i,-!  ,  013E*0'+ ) 

C0VALlEHALL2<f  ,EP      ,-1  .OlSE+O't.-l.OlSE+O** ) 


13,      1^, 


1) 


PATCH!  l'iWALL2<*  ,WWALL  ,  10,  10,  1,  8, 
C0VAL!WWALL2<^  ,V1  ,-1 .  013E  +  0<.,  0.000E*00> 
C0VAL!^T'ULL2^    ,W1      ,-1.013E+0^,    O.OOOE+OOJ 


13,      1<»,        1,        1) 
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C0VAL1WWALL2<»    ,KE 
C0VAL(WWALL2^    ,EP 


,-1.013E*0'»,-1.013E*0<i) 
,-1.013E*0«t,-1.013E+0^) 


PATCH(EWALL2S  ,EWALL    >        9,        9,        1,        8, 

C0VAL«EWALL2S  ,V1   ,-1.013E+0<*,  0.000E*00) 

COVAL(EHALL25  ,W1   ,-1.013E*0^,  O.OOOE+00) 

C0VAL(EHALL25  ,KE   ,-1 .  OlSE+O^.-l.OlSE+O'*  ) 

C0VAL(EWALL2S  ,EP   ,-1 .013E*0<»,-1.013E»0'»  ) 


PATCH (W-JALL2S  ,WWALL 

C0VAL(WWALL2S  ,V1   , 

C0VAL(HWALL2S  ,m   , 

C0VAL(WWALL2S  ,KE   , 

C0VAL(W>iALL2S  ,EP 


10,  10,  1,  8, 
,-I.013E*0<*,  O.OOOE+OO) 
,-1.015E+0<»,  O.COOE+00) 
,-1.012E*0'*,-1.013E*0<*) 
, -1 .  013E  ♦0'+ , -1 .  013E  ♦0«»  ) 


PATCH(EWALL26  ,EWALL  ,    5,    5,    1,    8, 

COVAL(EWALL26  ,V1   ,-1.013E+0<+,  0,000E*00) 

C0VAL(EWALL26  ,H1   ,-1  .OlSE+O'*,  O.OOOEtOO) 

COVAL(EWALL26  ,KE   ,-1.013E  +  0<+,-1.013E+0'» ) 

C0VAHEHALL26  >EP      ,-1 .013E+0<t  ,-1.013E+0't ) 


18,      19,        1, 


18,      19,        1, 


18,      39,        1, 


1) 


1) 


1) 


PATCH(WWALL26  ,WHALL 


6,        6,        1,        8,      18,      39, 


1, 


1) 


C0VAL(HWALL26  ,V1 
COVALIWWALL26  ,H1 
C0VAL(HWALL26  ,KE 
C0VAL(WWALL26  ,EP 
XCYCLE  =  F 
KELIN    =      0 

ifr^ic  Jt  w^iMtirTt  w  w  WTtitw^M*  ^^^^^^^M^^^  itHflMnvlviftMHMAIIK  fWUHtK  11  Jl  H  M  M  JWWUHHWm  M  N  M  H 

GROUP  l<i.  DOWNSTREAM  PRESSURE  FOR  PARAB 


,-1.013E+0^,  O.OOOE+00) 
,-1.013E*0<*,  O.OOOE+00) 
,-1.013E+0<f,-1.013E  +  0<+) 
,-1.013E*0<;,-1.013E+0'iJ 


GROUP    IS.    TERMINATE   SWEEPS 


w  ikiH^ik  ^  w  'M  R  wifr^  K  M  M  A  K  ff  ^Hr^r^f  it 


LSHEEP      = 

180 

jISWCl   = 

1 

LITHYD   = 

1 

JLITFLX  = 

1  JLITC    = 

1  )ITHC1 

RESREFIPI 

)  = 

l.OOOE-06 

JRESREF(U1   )  = 

1 

OOOE-06 

RESREFIVI 

)  = 

l.OOOE-06 

}RESREF(H1   )  = 

1 

OOOE-06 

RESREFIKE 

)  = 

l.OOOE-06 

>RESREF(EP      J    = 

1 

.OOOE-06 

RESREF(H1 

)  = 

l.OOOE-06 

GROUP    16.    TERi'UNATE    ITERATIONS 
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LITER  (PI 
LITER  I  VI 
LITER  (KE 
LITER  (HI 
ENDIT  (PI 
ENDIT  (VI 
ENDIT  (KE 
ENDIT  (HI 


20  iLITER  (Ul   )  =  20 

20  )LITER  (Wl   )  =  20 

20  jLITER  (EP   )  *  20 

20 

l.OOOE-06  }ENOIT  (Ul  )  s 

l.OOOE-06  jENDIT  (Wl  )  = 

l.OOOE-06  iENDIT  (EP  )  = 


l.OOOE-06 
l.OOOE-06 
l.OOOE-06 


=   l.OOOE-06 


fcj  W  ^J  W  W  V.  ^W^J^*JW^  w«->^>»J  'r    •,     W  \J   ^J^M 

H  M  n  M  ^W  W  .W  W  W  W  ^T^  H  K  M  ^^^^^  ^^^ 

GROUP  17.  RELAXATION 

RELAX(P1  ,LINRLX,  3.000E-01) 

RELAX(U1  ,FALSDT, 

RELAXCVl  ,FALSDT, 

RELAX(W1  ,FALSDT, 

RELAXIKE  ,FALSDT, 

RELAXIEP  ,FALSDT, 

RELAX(H1  ,FALSDT, 


RELAX(VIST,FALSDT, 
RELAX! RH01,FALSDT, 
RELAX(TMP1,FALSDT, 

GROUP  18.  LIMITS 


<».321E-05) 
<».321E-03) 
4.221E-03) 
321E-03) 
321E-03) 
080E-02) 
OOOEtOO) 
OOOE^OO) 
OOOEtOO) 


VARMAXfPl 
VARMAX(U1 
VARMAXl  VI 
VARMAX( Nl 
VARMAX(KE 
VARMAXl EP 
VARMAXl HI 
VARMAXl HPOR) 
VARMAXl NPOR ) 
VARMAXl EPOR) 
VARMAXl VPOR) 
VARMAXl VIST) 
VARMAXl RHOl) 
VARMAXl TMPl) 


3.500E+0S  jVARMINIPl 
l.OOOE+03  )VARMIN1U1 
1.000E+03  iVARMINlVl 
1.500E+03  JVARMINIHI 
l.OOOE+10  >VARMIN1K£ 
l.OOOE+10  JVARMINIEP 


)  =-3.S00E*05 
)  S-3.000E+02 
)  =-3.000E*02 
J  =-1.000E*02 
]  =  l.OOOE-10 
)  =  l.OOOE-10 


)  =  l.OOOE+03 


=  1. 


=  1.9S2E+06  iVARMINlHl 

=  1.000E*10  JVARMINIHPOR)  =-1.000E*10 

=  l.OOOE+10  iVARMINlNPOR)  =-1.000E+10 

=  l.OOOE+10  jVARMINlEPOR)  =-1.000E+10 

OOOE+10  JVARMINIVPOR)  =-1.000E*10 

l.OOOE+10  iVARMINlVlST)  =  l.OOOE-20 

S.OOOE+00  jVARMINl RHOl )  =  5. CODE -02 

1.9'»<*E+03  jVARMINl  TMPl )  =  2.500E+02 

GROUP  19.  EARTH  CALLS  TO  GROUTiD  STATION 
USEGRD   =    T   iUSEGRX  =    T 
NAMGRD   =NONE 
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GEW    = 

T 

GROUP  20. 

PRELIMINARY  PRINTOUT 

ECHO    = 

T 

GROUP  21. 

PRINT-OUT  OF  VARIABLES 

INIFLD   = 

F   )SUBWGR  = 

F 

OUTPUTIPI 

,Y,N,N,Y,Y,Y) 

OUTPUTCUl 

,Y,N,N,Y,Y,Y) 

OUTPUT( VI 

,Y,N,N,Y,Y,Y) 

GUTPUTIWl 

,Y,N,N,Y,Y,Y) 

DUTPUT(KE 

,Y,N,N,Y,Y,Y) 

OUTPUKEP 

,Y,N,N,Y,Y,Y) 

OUTPUTIHI 

,Y,N,N,Y,Y,Y) 

OUTPUT( HPOR,N,N,N,N,N,N) 

OUTPUT(NP0R,N,N,N,N,N,N) 

OUTPUT(EPOR,N,N,N,N,N,N) 

OUTPUT! V?OR,N,N,N,N,N,N) 

OUTPUT! VIST, Y,N,N,N,N,N) 

OUTPUT! RH01,y,N,N,N,N,N) 

OUTPUT! TMP1,Y,N,N,N,N,N) 

GROUP  22.  MONITOR  PRINT-OUT 

IXMON    = 

1  ilYKON   = 

<» 

lIZMON   = 

16 

NPRMON   = 

10000  >NPRMNT  =  10000 

)TSTSHP  = 

1 

HIGHLO   = 

F 

MMICII 

GROUP  23. 

FIELD  PRINT-OUT  i 

PLOT 

■  CONTROL 

NPRINT   = 

10000  iNUMCLS  = 

10 

NXPRIN   = 

1  ilXPRF   = 

1 

jIXPRL   = 

10 

NYPRIN   = 

1  JIYPRF   = 

1 

jIYPRL   = 

12 

N2PRIN   = 

1  JIZPRF   = 

1 

jIZPRL   = 

<*0 

X2PR     = 

F   }Y2PR    = 

F 

IPLTF    = 

1  ilPLTL   = 

180 

jNPLT    = 

1 

ITABL    = 

5  iIPROF   = 

1 

ABSIZ    = 

l.OOOE+00  J0RSI2 

=  1. 

OOOEtOO 

NT2PRF   = 

1  )NCOLPF  = 

50 

ICHR     = 

2  jNCOLCO  = 

tfS 

jNROWCO  = 

20 

NO  PATCHES  YET  USED  FOR  THIS  GROUP 
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GROUP  2<».  DUMPS  FOR  RESTARTS 
SAVE    c    T   jAUTOPS  =    F 
NSAVE    =DF09 
LUSAVE   =     9 


JNOWIPE   = 


«*«  GRID-GEOMETRY  INFORMATION  *** 

X- 

-COORDINATES 

OF  THE  CELL 

CENTRES 

1.S50E-01 

S.837E-01 

1.131E+00 

1 

.679E+00 

2. 

,226EtO0 

3.500E+00 

5.S00E+00 

7.S00E+OO 

9 

,200E+00 

1. 

.060E+01 

Y- 

-COORDINATES 

OF  THE  CELL 

CENTRES 

2.810E-01 

8.^50E-01 

l.<f05E*00 

1 

.843E+00 

2 

,lS7E+00 

2.59SE+00 

S.1S7E+00 

3.719E+00 

4 

.37SE*00 

5. 

,12SE*00 

S.87SE*00 

6.625E+00 

Z- 

-COORDINATES 

OF  THE  CELL 

CENTRES 

S.OOOE-01 

1.850E+00 

3.700E+00 

S 

.SSOE+OO 

7, 

,900Et00 

1.090E+01 

1.51SE+01 

1.<»<«0E+01 

1. 

.S75E+01 

1, 

.710E+01 

1.810E+01 

1.9'*Sc+01 

2.080E+01 

2 

.180E»01 

2. 

,258E*0I 

2.il5E+01 

2.372E+01 

2.417E+01 

2. 

,<V52E*01 

2. 

,620E*01 

2.920Et01 

3.220E+01 

3.S20Et01 

3. 

.820Et01 

<♦. 

.120E+01 

'♦.420E+01 

'».670E*01 

4.807Et01 

<t, 

.8a2E*01 

<t, 

.957E*01 

S.022E+01 

5.099E+01 

S.1S6E+01 

5, 

,21'*E*01 

S. 

,271E*01 

S.325E+01 

S.37SE*01 

5.<»2SE*01 

5, 

,'»75E+01 

5. 

bSOE+Ol 

INTEGRATION  OF   EQUATIONS  BEGINS   — ■ 


TIME  STEP  =    1 

SWEEP  = 

2 

TOTAL  RESIDUAL/! 

.OOOE-06) 

FOR 

PI 

IS 

8, 

,007E+0S 

TOTAL  RESIDUAL/( 

,000£-06) 

FOR 

Ul 

IS 

9, 

,092E+06 

TOTAL  RESIDUAL/( 

,000E-06) 

FOR 

VI 

IS 

1, 

.282E+07 

TOTAL  RESIDUAL/( 

.OOOE-06) 

FOR 

Wl 

IS 

2 

,376E*07 

TOTAL  RESIDUAL/( 

.OOOE-06) 

FOR 

KE 

IS 

1. 

.210E+09 

TOTAL  RESIDUAL/( 

.OOOE-06) 

FOR 

EP 

IS 

1 

.<»9eE  +  ll 

TOTAL  RESIDUAL/I 

.OOOE-06) 

FOR 

HI 

IS 

S 

.642E+10 

TIME  STEP  =    1 

SWEEP  = 

3 

TOTAL  RESIDUAL/I 

.OOOE-06) 

FOR 

PI 

IS 

3 

.733E+05 

TOTAL  RESIDUAL/I 

.OOOE-06) 

FOR 

Ul 

IS 

1 

.050E+07 

TOTAL  RESIDUAL/( 

.OOOE-06) 

FOR 

VI 

IS 

1 

.<»7'tE+07 

TOTAL  RESIDUAL/1 

.OOOE-06) 

FOR 

Wl 

IS 

2 

.037E+07 

TOTAL  RESIDUAL/( 

,0C0E-06) 

FOR 

KE 

IS 

1 

.'♦2<+Et09 
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TOTAL  RESIDUAL/( 
TOTAL  RESIDUAL/( 
TIME  STEP  =  1 
TOTAL  RESIDUAL/I 
TOTAL  RESIDUAL/( 
TOTAL  RESIOUAL/( 
TOTAL  RESIDUAL/( 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/( 
TOTAL  RESIDUAL/( 
TIME  STEP  =  1 
TOTAL  RESIDUAL/( 
TOTAL  RESIDUAL/( 
TOTAL  RESIDUAL/( 
TOTAL  RESIOUAL/( 
TOTAL  RESIDUAL/( 
TOTAL  RESIDUAL/( 
TOTAL  RESIDUAL/( 
TIME  STEP  =  1 
TOTAL  RESIDUAL/( 
TOTAL  RESIDUAL/t 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 
TIME  STEP  =  1 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 
TIME  STEP  =  1 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 
TOTAL  RESIDUAL/l 


OOOE-06) 

FOR 

EP 

IS 

2 

.3<*1E*10 

OOOE-06) 
SWEEP  = 

FOR 
17<* 

HI 

IS 

2 

.620E+10 

OOOE-06) 

FOR 

PI 

IS 

3 

231E+04 

OOOE-06) 

FOR 

Ul 

IS 

6 

021E+06 

OOOE-06) 

FOR 

VI 

IS 

6 

189E+06 

OOOE-06) 

FOR 

HI 

IS 

8 

271E+06 

OOOE-06) 

FOR 

KE 

IS 

2 

926E*0a 

OOOE-06) 

FOR 

EP 

IS 

2 

^5£E*10 

OOOE-06) 
SWEEP  = 

FOR 
175 

HI 

IS 

2 

619E+10 

OOOE-06) 

FOR 

PI 

IS 

3 

208EtO<i 

OOOE-06  ) 

FOR 

Ul 

IS 

5 

996E*06 

OOOE-06) 

FOR 

VI 

IS 

6 

183Et06 

OOOE-06) 

FOR 

Wl 

IS 

8 

329E+06 

OOOE-06) 

FOR 

KE 

IS 

2 

91SE+08 

OOOE-06) 

FOR 

EP 

IS 

2 

395Etl0 

OOOE-06) 
SWEEP  = 

FOR 
176 

HI 

IS 

2 

628E^10 

OOOE-06) 

FOR 

PI 

IS 

3 

193E*0<* 

OOOE-06) 

FOR 

Ul 

IS 

6 

071Et06 

OOOE-06) 

FOR 

VI 

IS 

6 

12SE*06 

OOOE-06) 

FOR 

HI 

IS 

8 

2S9E+06 

OOOE-06) 

FOR 

KE 

IS 

2 

909E+0a 

OOOE-06) 

FOR 

EP 

IS 

2 

-♦SSE+IO 

OOOE-06) 
SWEEP  = 

FOR 
177 

HI 

IS 

2 

628E+10 

OOOE-06) 

FOR 

PI 

IS 

3 

235E*0'+ 

OOOE-06) 

FOR 

Ul 

IS 

6 

038E+06 

OOOE-06) 

FOR 

VI 

IS 

6 

128E+06 

OOOE-06) 

FOR 

Wl 

IS 

8 

291E+06 

OOOE-06) 

FOR 

KE 

IS 

2 

9S'»E*08 

OOOE-06) 

FOR 

EP 

IS 

2 

'♦66E  +  10 

OOOE-06) 
SWEEP  = 

FOR 
178 

HI 

IS 

2 

622E+10 

OOOE-06) 

FOR 

PI 

IS 

3 

272E+04 

OOOE-06) 

FOR 

Ul 

IS 

6 

09SE+06 

OOCE-06) 

FOR 

VI 

IS 

6 

115E+06 

OOOE-06) 

FOR 

HI 

IS 

8 

309Et06 

OOOE-06) 

FOR 

KE 

IS 

3 

OlOE+08 
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TOTAL  RESIDUAL/( 
TOTAL  RESIDUAL/I 
TOTAL  RESIDUAL/( 
TOTAL  RESIOUAL/t 
TOTAL  RESIDUAL/I 


l.OOOE-06)  FOR  VI    IS   6. 
l.OOOE-06)  FOR  Wl    IS   8, 
l.OOOE-06)  FOR  KE    IS   2. 
l.OOOE-06)  FOR  E?        IS   2, 
l.OOOE-06)  FOR  HI    IS   2. 
WHOLE-FIELD  RESIDUALS  BEFORE  SOLUTIONS 
WHOLE-FIELD  SUM  OF  ABS( VOL. FLOW  RESIDUALS)' 
WHOLE-FIELD  SUM  OF  ABS( RESIDUALS )  OF  Ul    ' 
WHOLE-FIELD  SUM  OF  ABS( RESIDUALS )  OF  VI 
WKOLE-FIELD  SUM  OF  ABSt RESIDUALS )  OF  Wl 
WHOLE-FIELD  SUM  OF  ABS( RESIDUALS )  OF  KE 
WHOLE-FIELD  SUM  OF  ABSI RESIDUALS )  OF  EP 
WHOLE-FIELD  SUM  OF  ABS(  RESIDUALS )  OF  HI    ■■ 
•   SUMS  HAVE  BEEN  DIVIDED  BY  RESREF(NAME) 


1<*4E«06 
IZIE*0(, 
930E+08 
^30E+10 
630Etl0 


S.185EtO<i 
6.088Et06 
6.1<^E«06 
8.323E*06 
2.950E+0a 
2.<t50E  +  10 
2.630EtlO 


NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOUI^CE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOURCE 
NET  SOUr?C£ 
NET  SOURCE 


OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 
OF  Ul 


AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 
AT  PATCH 


NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED: 
NAMED; 


NWALL3 

SWALL3 

NWALL<» 

SHALL<^ 

NWALLS 

SWALLS 

NWALL6 

SHALL6 

NWALL7 

SWALL7 

NWALLS 

SWALLS 

SWALL9 

HWALLIO 

LWALLIO 

HWALLll 

LWALLll 

HWALL12 

LWALL12 

HWALL13 

LWALL13 

HWALLll 

LWALL1<+ 

HWALLIS 

LWALL15 


=  -<j.696E*00 
=  1.728E-01 
S-1.275E-02 
S-3.272E-01 
=-2,179£-02 
=-1.7S7E-01 
=-2.196E-02 
=-2.1^';E-01 
S-1.102E-02 
=-1.9bOE-01 
=-6.52^E-03 
S-1.813E-01 
=  2.277E+01 
»  O.OOOE+00 
=  8.S69E-01 
=  1.661E+00 
=-2.399E-02 
S-9.981E-02 
=-6.37'^E-0'V 
s-l.OOOE+00 
=-1.89SE-02 
S-6.707E-01 
=-2.^7=£-02 
=-1.2b^E-01 
=-1.210E-02 
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NET 

SOURCE 

OF 

Ul 

AT 

PATCH 

NAMED: 

HWALL16 

S' 

-3 

.<»65E-01 

NET 

SOURCE 

OF 

Ul 

AT 

PATCH 

rUHED: 

LWALL16 

S' 

-1, 

.209E-02 

NET 

SOURCE 

OF 

Ul 

AT 

PATCH 

NAMED: 

HWALL17 

s- 

-3 

,<t2^E-01 

NET 

SOURCE 

OF 

Ul 

AT 

PATCH 

NAMED: 

LWALL17 

S' 

-5, 

.672E-03 

NET 

SOURCE 

OF 

Ul 

AT 

PATCH 

NAMED: 

HKALL18 

r. 

-3, 

.36<tE-01 

NET 

SOURCE 

OF 

Ul 

AT 

PATCH 

NAMED: 

LWALL18 

= 

8, 

,109E-02 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

HMALLIO 

: 

0, 

.OOOE+00 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

LWALLIO 

S' 

-1. 

,138E+00 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

HHALLll 

= 

2, 

.250E+00 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

LWALLll 

= 

9, 

.259E-03 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

HHALL12 

= 

2, 

.302E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

LWALL12 

s 

<►, 

.28'^E-O't 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

HWALL12 

=  . 

-9, 

,006E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

LWALL13 

s 

1. 

,062E-03 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

HWALLl'* 

s- 

■1, 

.16HE+00 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

LNALLl'i 

s 

1. 

,3<t2E-0<+ 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

HWALLIS 

=  • 

-2, 

,973E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

LWALLIS 

s. 

-<♦, 

.722E-07 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

HWALL16 

s. 

-8. 

.^32E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

LWALL16 

=  ■ 

-6, 

.630E-07 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

HHALL17 

s- 

-9. 

,907E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

LHALL17 

=■ 

-6, 

,93^E-07 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

HWALLIS 

r- 

-1. 

.293E+00 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

LWALLIS 

s- 

-1, 

.613E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

EWALL19 

s 

8, 

.668E-01 

NET 

SOURCE 

Or 

VI 

AT 

PATCH 

NAMED: 

WWALL19 

s 

3. 

.758E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

EWALL20 

s- 

-8. 

.752E-03 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

WWALL20 

= 

3. 

.321E-02 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

EWALL21 

= 

.2S4E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAJ1ED: 

WWALL21 

s 

,022E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

EWALL22 

s 

.8'+7E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

KWALL22 

: 

,^80E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

EWALL23 

= 

,075E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

WWALL25 

= 

.80SE-02 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

EWALL2<f 

= 

.69<+E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

NWALL2<+ 

: 

,0^5E-02 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

EWALL2S 

s 

,168E-01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

WWALL25 

= 

,S18E-02 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

EWALL26 

=  ■ 

-2 

.0'+7E+01 

NET 

SOURCE 

OF 

VI 

AT 

PATCH 

NAMED: 

WWALL26 

s- 

-3, 

.99SE-01 
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NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

JETOUT 

s 

^ 

.176E*0<* 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

NHALL5 

S' 

-1. 

.318E»02 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

SHALL3 

S' 

-1 

.018E+00 

NET 

SOURCE 

OF 

Wl 

AT 

PATCH 

NAMED: 

NWALL<i 

s 

3 

.301E-03 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

SHALL^ 

S' 

-1 

.247E-08 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

MiALLS 

s 

1 

.930E-03 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

SHALL5 

S' 

-8 

.121E-09 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

NHALL6 

S' 

-3 

.027E-03 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

SHALL6 

S' 

-1 

.<*50E-08 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

NWALL7 

S' 

-8 

.7'«3E-03 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

SHALL7 

=• 

-1, 

.S02E-08 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

NHALL8 

=' 

-1, 

,9SlE-02 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

SHALL8 

=  • 

-1. 

.S^OE-03 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

SHALL9 

S' 

-8, 

,925E+01 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

EHALH9 

s- 

-2. 

,025E^00 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

HHALL19 

s 

6, 

.S76E-02 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

EHALL20 

s 

6. 

.667E-02 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

HVIALLZO 

=■ 

-1, 

.278E-01 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

EHALL21 

= 

2. 

.<»22E-01 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

HHALL21 

:• 

-1. 

.12^E-01 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

EHALLZ2 

*• 

-2, 

.738E-01 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

HHALL22 

S' 

-3. 

.260E-01 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

EHALL23 

s 

1. 

.531E-02 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

KHALL25 

s< 

-1, 

.3S0E-01 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

EHALL2«* 

s 

9, 

,^l<*E-02 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

HHALL2<;^ 

s> 

-<t. 

,S38E-02 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

EHALL25 

r 

<*. 

.882E-0^ 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

HWALL25 

s 

1, 

,093E-01 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

EKALL26 

s< 

-2. 

,526E+01 

NET 

SOURCE 

OF 

HI 

AT 

PATCH 

NAMED: 

HHALL26 

S' 

-1. 

.1<»0E-01 

NET 

SOURCE 

OF 

Rl 

AT 

PATCH 

NAMED: 

JETIN 

z- 

-3, 

,8<»9E*01 

NET 

SOURCE 

OF 

Rl 

AT 

PATCH 

NAMED: 

JETOUT 

s 

3, 

,6^9E+01 

NET 

SOURCE 

OF 

Rl 

AT 

PATCH 

NAMED: 

OUTSIDE 

: 

<*. 

,1S5E*02 

NET 

SOURCE 

OF 

Rl 

AT 

PATCH 

NAMED: 

TOPO'JT 

=  • 

-<*. 

.OlOE+02 

NET 

SOURCE 

OF 

Rl 

AT 

PATCH 

NAMED: 

TOPP 

= 

2, 

,728E+01 

NET 

SOURCE 

OF 

Rl 

AT 

PATCH 

NAMED: 

FRONT 

= 

1. 

.8S9E+01 

NET 

SOURCE 

OF 

Rl 

AT 

PATCH 

NAMED: 

ENDEND 

S' 

-6, 

.022E+01 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

RESOURCE 

= 

6. 

.lOOE+05 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

JETOUT 

= 

5. 

.663E«02 

NET 

SOURCE 

OF 

K£ 

AT 

PATCH 

NAMED: 

NWALL3 

r. 

-1. 

,089E+10 
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NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

SWALL3 

*■ 

-1 

.308E«10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

NWALL<» 

= 

3 

.850E«09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

SWALL<» 

=  • 

-8 

.IIIE^IO 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NA/1ED: 

NUALLS 

= 

<i. 

.497E*09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

SWALLS 

=• 

-6, 

,297E<-10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

NKALL6 

: 

<i. 

.516Et09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

SWALL6 

S' 

-1, 

.9«E  +  11 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

NWALL7 

1 

.865E*09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

SWALL7 

-2 

.<»37E  +  11 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

Ni'iALLS 

6, 

.060E+09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

SHALL8 

-2 

.84SE+11 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

SWALL9 

-6, 

,320E*10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

HWALLIO 

0 

.OO0E*00 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

LWALLIO 

-2. 

.OOSE+10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

HWALLIl 

= 

8, 

.3<t4EflO 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

LWALLll 

=  • 

-3 

,S67E*Q8 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

HKALL12 

= 

1. 

,310E*10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

LWALL12 

=  • 

-1. 

,620E+07 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

HWALLIJ 

s 

7, 

,1S5E+10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

LWALL15 

s- 

-8, 

,321E+09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

HWALLl<t 

s 

9. 

,199Et-10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

LWALLl<i 

=  • 

-4, 

,643E+09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

HWALL15 

s 

6. 

,299E+10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

LHALLIS 

:• 

-<♦. 

.610E+09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

HWALL16 

s 

1. 

.840E+11 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

LWALL16 

s- 

-2, 

.695E+09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

HWALL17 

= 

2. 

.26SE+11 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

LWALL17 

s< 

-7, 

,609E+09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

HWALLia 

s 

2. 

,606E+11 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

LWALL18 

s 

6. 

,3S8E+05 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

EWALL19 

=  ■ 

-<♦. 

.326E+09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

HWALL19 

= 

3. 

,  29'4E*05 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

EWALL20 

s 

6. 

.926E+09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

KHALL20 

=■ 

-8, 

,3'+9£*08 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

EWALL21 

=  • 

-2, 

,170E+09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

WWALL21 

=  • 

-3, 

,  718E*08 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

EHALL22 

s- 

-2, 

.S67E+09 

NET 

SOURCE 

Or 

KE 

AT 

PATCH 

NAMED: 

WHALL22 

s< 

-2, 

,501E*09 

NET 

SOURCE 

Or 

KE 

AT 

PATCH 

NAMED: 

EHALL23 

=  ' 

-<♦, 

,046E+09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

HKALL22 

=  • 

-9. 

,862E^0d 
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NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

EWALL24 

=  -<♦, 

,0<»5E*09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

WHALL2<» 

s-2 

.089E+08 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

EHALL2S 

=-1, 

.601E*10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

WWALL2S 

--Z. 

,325E+09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED: 

EKALL26 

=  8. 
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AJPPEa^DDIX    C 


AIRFLOW  CHARACTERIZATION  SHEET 


TEST  PROCRAMi      Tesclne  of   Navy  Aircraft    In  Arched   Structure   Hush   House 
tOCATlOH.  K^lly   A.F.B 


OATEi       9    November    1982 
AIRCRAFT:        F-AN 


ewQIHEi      J-79-CE-8 
EVENT{  5 


CELL    DEPRESSION    (l«.    WB.K       1-7 


CONDITIONS!      Port   Engine   -  A/B   Power.    Starboard   Engine  -   Idle   Power 

AUOMENTER   INLET 


:}0  ■•;:<.'0-c;::-:- 


a 

;,,^y,/,,-./^i.\y,-,; 


i 


Limited  D«ca 
Recorded 


■  ' '•■>  *'     *■■■  ■*-■■■■■        ■  • »»    '■%■,*-"*'_  \*       .  I     ''--'i    ■■■'■■^ 


^ 


12.7   f.i 


MUSH    HOUSE    PLAN    VIEW 


126 
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AIRFLOW  r.HARACTERIZATION  SHEET 


;"^"/^;i°.°"^''N"val    Air    Stat^on^    .Tack.sonvi  1  1  p      Flnfil^ 
DATE:   _. 


AIRCRAFT: 


N/A 
N/A 


gMciMEi     J-79-GE-8 
N/A 


Cell  depression  (i*.  »■•>-     ^/a 


EVENT:   , . 

Sinqle  Engine  with  A/B  (g  rpnt.Prlme 

conditions:     "^  "^  ^  ^ ' 

auouenterinlet 


o 


Ik 
o 


z 
o 


K 

o 


TVPtCAt 
.|.0W>EC0CIJf 


/^mPLOW'BEOlOlifi 


HUSH    HOU 


8E    PLAN    VIEW  1 
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Ar>PENI>IX    D 


EXHAUST  FLOW  DATA  SHEET 


TEST  PROCRAM:       Testing  of  Wavy  Aircraft   in  Arched   Structure   Hush  House 
LOCATION:        Kelly  A.F.B.  ,    San  (Lotoplo,   Texas 


OAXE:       9  November  1962 
AIRCRAFT:        F-4-N 


EMGINE;    J-79-CE-8 
EVENT:   S 


CONDITIONS:         Port   Engine  -  A/B,    Starboard   Engine  -   Idle   Power 


AUCMENTER    TUBE 


•  (•rboard 


V=VELOCiTY    (feet/sec.) 


V=    209.5 

• 

V  = 

• 

T= 

271.6 

A54 

V=  239.7 

• 

T=454 

V=    20i.9 

• 

T=    <i5<. 

v= 

• 
T= 

29A.8 
434 

V=    372.9 

• 

T=  454 

V=    156.5 
T=  454 


V=    150.1 
T=    454 


T=TcMPERATURE  (degrees  F) 


TOTAL   MASS   FLOW   (Ibm. /••€.)=     2067 
AUCMENTATIOM    RATIOS    8.3 


Vavg=  ZZl.Z 
Tavg=  45^ 


part 
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\ 


EXHAUST  FLOW  DATA  SHEET 


TEST  PROGRAM:        TheSJS 


LOCATION:    Naval   Air   Station,    Jacksonville.    FlnriH;:) 

DATE:   .    N/A 


"nTa" 


AIRCRAFT:  

ewGiME:   J-79-GE-8 

EVENT:  ^^^ 


CONDITIONS:    Single   Engine   with  A/B    P    centerlinp 


AUGMENTER    TUBE 


•larboard 


V=V£LOCITY   (feet/sec.) 

v=236  v=305.9  v= 


V=167.3 
T=  395.6 


T=426.2 

V=318.3 
T=453.2 


V=233 

T=386.6 


T= TEMPERATURE  (degrees  F) 


FLOW  (lbm./»«e.)s  237  3 
AUGMENTATION  RATIO"  10.6 


Vavg=    235.5 
Tavg=    401.6 


pert 
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APPENDIX  E 
DiATA  TRANSFER  BY  TAPE 


The  IBM  3033  main-frame  presently  does  not  have  "PHOTON" 
graphics  capabilities.  In  order  to  produce  various  visual 
plots  that  enhance  the  interpretation  of  tabulated  data,  the 
"DF09"  and  "DF12"  files  must  be  properly  formatted  and  then 
transferred  to  tape  (for  portability)  allowing  use  of  "PHOTON" 
on  another  computer  which  is  "PHOTON"  capable. 

The  computer  center  will  provide  a  non-labelled  tape  for 
data  transfer.  The  operator,  when  queried  by  the  programs 
below,  will  mount  the  tape  allowing  the  "DF09"  and  "DF12" 
files  to  be  transferred  to  the  tape. 

To  prepare  the  files  for  transfer  to  tape  a  copy  is  made 
to  ensure  that  the  original  data  file  is  safe.  Program  "FORM09 
JCL"  properly  formats  and  then  writes  the  "DF09"  file  to  a  new 
file  called  "DF09C".  Program  "F0RM12  JCL"  properly  formats  and 
then  writes  the  "DF12"  file  to  a  new  file  called  "DF12C".  The 
data  sets  "DF09C"  and  "DF12C"  which  are  produced  using  MVSHELP 
must  be  at  least  as  large  as  the  original  "DF09"  and  "DF12" 
data  sets  and  created  before  running  "F0RM12  JCL"  or  "FORM09 
JCL". 

With  the  two  files  now  properly  formatted  the  file  "DSTP 
JCL"  is  used  to  transfer  the  designated  file  to  tape.  In  the 
following  example  the  "DF09C"  file  is  transferred  as  the 
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second  file  to  the  tape  called  "NICOLA".  File  "TSCAN  JCL"  then 
scans  the  tape  to  ensure  the  program  was  actually  transferred 
to  the  tape. 

Two  other  files,  "TAP  JCL"  and  "DUMP  JCL",  can  dump  a 
file  from  tape  to  the  virtual  machine  and  return  the  first 
block  of  the  tape  file  respectively.  With  the  files  on  tape 
and  checked  the  tape  can  be  dismounted  by  the  operator  at  the 
computer  center  and  transported  to  a  "PHOTON"  capable  system 
in  the  "formatted"  form. 

The  various  files  referenced  are  now  listed. 


FILE:  DSTP     JCL      Al 


//DSTP  JOB  (IS^tl, 9999), -ERIC  •,CLASS=E 

//   EXEC  PGM=IEBGENER 

//SYSPRINT   DD   SYSOLrr=A 

//SYSIN  DD  DUMMY 

//SYSUTl  DD  DISP=SHR,OSN=MSS. 515^1. PHOENICS.DF09C 

//SYSUT2  DD  UNIT=3^00-6,VOL=SER=NICOLA,DISP=(NEN,PASS), 

//        LABEL=(2,NL,,0UT), 

//       DCB=(RECFM=FB,LRECL=80,BLKSIZE=800,DEN=2,OPTCD=Q) 
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FILE:  TAP      JCL      Al 


//TAP  JOB  (1541,9999), 'TAP  JCL',CLASS=E 

//»MAIN  CAROS=(SO) 

//   EXEC  PGM=IEBGENER 

//SYSPRirfT  DD  SYSOUT=A 

//SYSIN  DD  DUMMY 

//SYSUTl  DD  UNIT=3<:*00-6,VOL=SER=NICOLA,DISP  =  (OLD,PASS), 

//       LABEL=(2,NL,,IN), 

//      DCB=(RECFM=FB,LRECL=80,BLKSIZE=800,OPTCD=Q) 

//SYSUT2  DD  SYSOUT=B,DCB=BLKSIZE=80 

// 


FILE:  DUMP     JCL      Al 


//DUMP  JOB  (1541,9999),' DUMP ',CLASS=E 

//  EXEC  TAPE,VOLIN=NICOLA 

//SYSIN  DD  » 

DMPRECC 10,272,1 ),SKPFIL( 10,1 ),XECLrrE(  999) 

/* 

// 


FILE:  TSCAN     JCL       Al 


//6SCAN  JOB  (1541, 9999),' NICOLA ',CLASS=E 

//   EXEC  TSCAN, VOLIN=NICOLA 

// 
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FILE:  FORM09   JCL      Al 


//FORM09  JOB  (15<>1, 9999), 'CHAM ',CLASS=B 

//«MAIN  LINES=(20) 

//  EXEC  FORTVCLG 

//FORT.SYSIN  DD  » 

C     PROGRAM  F0RM09 

C 

C  Increase  PARAMETER  settings  if  more  variables  or  a  larger  grid 

C  is  required. 

C 

PARAMETER( NPHIP=80 ) 
C  NB  NXYMAX  musi  be    >=   HZ   and  NX*NY 

PARAMETER!  h4XYMAX=20000  ) 

CHARACTER*"^  NAME!  NPHIP  )  ,MESS(10  ) 

CHARACTER*15  FILIN,FILOUT 

CHARACTER*!  ANS 

CHARACTER«10  ACCESS 

LOGICAL  LBUFl 10  ),STORE(NPHIP), DEBUG 

DIMENSION  RBUF( NXYMAX  ),IBUF (10) 
C 

100  FORMAT (IX, 71 10) 

101  F0RMAT(6(1PE13.6)) 

102  F0RMAT(1X,79L1) 
105  FORMAT (IX, 19A^) 

C 

C  Resci  DEBUG=.TRUE.  -to  gei   error  messages.. 

DEBUG=.TRUE. 
C 
C  Sei  Logical  Uniis... 

LUIN=9 

LU0UT=10 

LUBUG=6 
C 

C  Set  File  Names... 
C  Name  of  unformatted  restarf  file 
C     WRITE (6,*)  '  RESTART  FILENAME?' 
C     READ(5,'(A)' )  FILIN 
C  Name  of  formatted  restart  file 
C     HRITE(6,»*)  '  FORMATTED  FILENAME?' 
C     READ(5,'(A)' )  FILOUT 
C 

C  Direct  access  or  sequential? 

C     WRITE(6,*)  '  DIRECT  ACCESS  (PHIDA)  FILE?  Y/N' 
C     READ(S,'(A1)' )  ANS 
C     IF(ANS.EQ. 'Y'.OR.ANS.EQ. 'Y' )  THEN 
C  Version  l.<»  default 
C     ACCESS=' DIRECT' 
C     LENREC=256 
C     ELSE 
C  Version  l.S 

ACCESS= 'SEQUENTIAL' 
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FILE:  FORM09   JCL      Al 


LENREC=100 
C     ENDIF 
C     MRITE(6,»)  •  DEFAULT  LENREC?  (100  VI. 5,  2S6  VI. <f)  -  Y/N' 

C     READ(5,'(A)'  )  ANS 

C      IFIANS.NE. 'Y' .AND.ANS.NE. 'Y' )  THEN 

.C     WRITE (6,*)  '  LENREC?  (WORDS  NOT  BYTES)' 

C     READ(5,'(I)' )  LENREC 

C      ENDIF 

C  NBYTES  is  ihe  number  of  bytes  per  single  precision  word. 

NBYTES=^ 
C  Record  lengih  in  bytes. 

LREC=NBYTES«LENREC 
C 
C   OPEN  inpui  file: 

OPEN! UNIT=9,F0RM= 'UNFORMATTED' ,STATUS='OLD ' , 
1  ACCESS=ACCESS,IOSTAT=IERR,ERR=90 ) 
IF(IERR.NE.O)  GO  TO  90 

IF(ACCESS.EQ. 'DIRECT' )  THEN 

READILUIN,REC=1)    ( MESSd  ),I=1,10) 

ELSE 

READ(LUIN)  (MESSd), 1=1, 10) 

ENDIF 

IF( DEBUG)  WRITE! LUBUG,*)  MESS 

IFIACCESS.EQ. 'DIRECT'  )  THEN 

READ(LUIN,REC=2)  (  LBUF(  I  ),I=1,<*  ) 

ELSE 

READ(LUIN)  (LBUF(I),I=1,^) 

ENDIF 

IF( DEBUG)  WRITE! LUBUG,«)  (  LBUF( I ),1=1,4 ) 

IF( ACCESS. EQ. 'DIRECT' )  THEN 

READ(LUIN,REC=2)  NX, NY, NZ, NPHI ,(IBUF(I ), 1  =  1,6  ), LENREC 

ELSE 

READ(LUIN}  NX,NY,NZ,NPHI,(IBUF(I),I=1,6),LENREC 

ENDIF 

NXNY=NX*NY 

IF(NXNY.GT.NXYMAX)  GOTO  92 
IF(NPHI.GT.NPHIP)  GO  TO  93 
IF(DEBUG)  WRITE(LUBUG,«)  NX,NY,N2 
C 
C   OPEN  output  file:  FORMATTED,  SEQUENTIAL  and  of  FIXED  RECORD  LENGTH. 

OPEN( UNIT=10 ,FORM= ' FORMATTED ' ,STATUS= 'OLD ' , 
1  ACCESS= 'SEQUENTIAL ' ,I0STAT=IERR,ERR  =  91 ) 
IF(IERR.NE.O)  GO  TO  91 

WRITE! LUOUT, 103)  ( MESS( I  )  ,1  =  1 ,10  ) 

WRITE!  LUOUT,  102)  (  LBUF!  I  )  ,1  =  1  ,<> ) 

WRITE  I LUOUT , 100 )  NX ,NY ,NZ ,NPHI ,  I IBUF 1 1 ) ,1  =  1 ,6 ) , LENREC 

IFIACCESS.EQ. 'DIRECT'  )  THEN 

READ(LUIN,REC=4)  RINfJER 

ELSE 

READ!LUIN)  RINNER 
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FILE:  F0RM09   JCL      Al 


ENDIF 

WRITE! LUOUT, 101)  RINNER 
IF( ACCESS. EQ. 'DIRECT' )  THEN 
READ(LUIN,REC=5)  (NAME( I )>I=1,NPHI ) 
ELSE 

READ(LUIN)  (NAME(I),I=1,NPHI) 
ENDIF 
C  Set  C56...  rtames  for  variables  SI  io  NPHI  for  version  1.4,  "to  avoid 
C  control  characters  etc  I  (Unless  already  stored) 
IF(ACCESS.EQ. 'DIRECT' .AND. NPHI .GT. 50 )  THEN 
READ(LUIN,REC=10)  C STOREt I  ),I=1,NPHI ) 
DO  6  1=51, NPHI 
J=I-15 
IF( .NOT.STORE(I))  THEN 

IF(I.LT.IOO)  THEN 

WRITE(NAME(I),'(  "C",I2,"  ")')  J 

ELSE 

WRITE(NAME(I),'C  "C",I3)' )  J 

ENDIF 
ENDIF 
6  CONTINUE 
ENDIF 
WRITE (LUOUT, 103)  INAME( I ) ,I=1,NPHI ) 

IF(DEBUG)  WRITE (LUBUG,*)  (NAME! I ), 1=1, NPHI ) 
IF(ACCESS.EQ. 'DIRECT' )  THEN 
READILUIN,REC=6)  (  RBUF( I ),I=1,NX) 
ELSE 

READ(LUIN)  (RBUF(I),I=1,NX) 
ENDIF 

WRITE(LUOUT,101)    ( RBUFd  ),I  =  1,NX) 
IF(ACCESS.EQ. 'DIRECT' )    THEN 
READ(LUIN,REC=7)    ( RBUF( I ),I=1,NY J 
ELSE 

READ(LUIN)  (RBUF(I),I=1,NY) 
ENDIF 

WRITE (LUOUT, 101)  ( RBUF( I ),I=1,NY ) 
IF(ACCESS.EQ. 'DIRECT' )  THEN 
READ(LUIN,REC=8)  (  RBUFd  ), 1  =  1, N2) 
ELSE 

READ(LUIN)  (RBUF(I),I=1,N2) 
ENDIF 

WRITE! LUOUT, 101)  ( RBUFt I ) ,I=1,NZ ) 
IFIACCESS.EQ. 'DIRECT' )  THEN 
READILUIN,REC=9)  ( RBUF( I ),I=1,NZ) 
ELSE 

READ(LUIN)  (RBUFd), 1=1, NZ) 
ENDIF 

WRITE  (LUOUT,  101)  (  RBUFd  )  ,I=1,NZ) 
IF(ACCESS.EQ. 'DIRECT' )  THEN 
READ(LUIN,REC=10)  ( STORE d ), 1=1, NPHI ) 
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FILE:  FORM09   JCL      Al 


ELSE 

READ(LUIN)  (ST0RE(I),I=1,NPHI) 

ENDIF 

HRITE( LUOUT , 102 )  ( STORE ( I ) ,1  =  1 ,NPHI ) 

IF(DEBUG)  HRITE(LUBUG,102)  (STORE( I ),1=1,NPHI ) 

NVAR=0 

DO  S  I=1,NPHI 

IF(STORE(I))  NVAR=NVAR+1 
5  CONTINUE 

IREC=10 

DO  10  K=1,NZ 

DO  20  I=1,NVAR 

IFIACCESS.EQ. 'DIRECT' )  THEN 

IREC=IREC+1 

NVREC=NXNY/LENREC 

IFI MODI NXNY ,LENREC ) .NE .0 )  NVREC=NVREC*1 

00  30  IVREC=1,NVREC 

JF=(IVREC-1)*LENREC  ♦  1 

JL=IVREC»LENREC 

JL=MINO( JL,NXNY) 

READ(LUIN,REC=IREC)  ( RBUF( J  )  ,J=JF,JL  ) 
50  CONTINUE 

ELSE 

READ(LUIN)  (RBUF(J),J=1,NXNY) 

ENDIF 

WRITE! LUOUT, 101)  (  RBUF( J ), J=1,NXNY  ) 
20  CONTINUE 
10  CONTINUE 

IF( DEBUG)  HRITE(LUBUG,1900) 

1900  FORMAT! IX,'  NORMAL  STOP  IN  PROGRAM') 
GO  TO  999 

C  Error  irapping. . . 

90  IF( DEBUG)  WRITE!  LUBUG,1901) 

1901  FORMAT! IX,'  ERROR  IN  OPENING  INPUT  FILE') 
GO  TO  999 

91  IF! DEBUG)  WRITE! LUBUG, 1911 ) 

1911  FORMAT!  IX, 'ERROR  IN  OPENING  OUTPUT  FILE') 
GO  TO  999 

92  IFINXNY.GT.NXYMAX)  WRITE!  LUBUG, 1921 ) 
IF1N2.GT.NXYMAX)  WRITE! LUBUG, 1925  ) 

1921  FORMAT! IX, 'SETTING  OF  NX*NY  TOO  LARGE)  RESET  IN  FORM09  AND 

1,/, IX, 'RE-CREATE  THE  LOAD  MODULE') 
1923  FORMAT! IX, 'SETTING  OF  HZ   TOO  LARGE)  RESET  IN  FORM09  AND  ' 
1,/, IX, 'RE-CREATE  THE  LOAD  MODULE') 
GO  TO  999 
95      WRITE! LUBUG, 195) 
195  FORMAT! IX, 'SETTING  OF  NPHI  TOO  LARGE)  RESET  IN  FORM09  AND  ' 

1,/, IX, 'RE-CREATE  THE  LOAD  MODULE') 
999  CONTINUE 
C   CLOSE  files  ... 
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FILE:    FORM09        JCL  Al 


CLOSE(LUIN,STATUS='KEEP' ) 

CLOSECLUOUT,STATUS='KEEP'  ) 

STOP 

END 
//GO.FT09F001  DD  DSN=MSS.S15'tl .PHOENICS.DF09,DISP  =  (OLD  ) 
//60.FT10F001  DD  DSN=MSS.S1541.PHOENICS.DF09C,DISP=( OLD ) 
/« 
// 
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FILE:  F0RM12   JCL      Al 


//F0RM12  JOB  (1541, 9999), 'CHAM ',CLASS=B 

//«MAIN  LINES=(20) 

//  EXEC  FORTVCLG 

//FORT.SYSIN  DD  • 

C  File  name:  F0RM12  lasi  amended  15  October  1986 

C     PROGRAM  F0RM12 

C Program  io  rewrite  a  DF12  file  into  a  formatted  READCO  file. 

C   (JPE  S  JFD  ) 

C 

C  Reset  NXYIP  if  larger  grid  required. . 

PARAMETER( NXY1P=10000 ) 

DIMENSION  AREC( NXYIP) 
C  Reset  fixed  record  length  of  DF12  file  if  required. . . . 

CHARACTER*30  FILIN,FILOUT 

LOGICAL  DEBUG 

INTEGER  SECTOR 

DATA  LENREC/100/ 
C 

100  FORMAT! 315) 

101  FORMAT! 5(1PE13.6)) 
C 

C  SECTOR  represents  the  largest  record  length  that  can  be  OPENed. 

SECT0R=65536 
C  MBYTES  is  the  number  of  bytes  per  single  precision  word. 

NBYTES=<+ 

LREC=MIN0(NBYTES*LENREC,SECTOR-NBYTES) 

IBLK=LREC+NBYTES 
C 
C  Set  Logical  Units... 

LUIN=12 

LU0UT=13 

LUBUG=6 
C  Reset  DEBUG  if  error  messages  required. . 

DEBUG=.TRUE. 
C  Set  file  names 

C      FILIN='DF12  • 

C     FILOUT='DF12F0RM  * 

C 

C  OPEN  input  file:  DIRECT  ACCESS,  UNFORMATTED  and  of  FIXED 
C  RECORD  LENGTH. 
C 

OPEN( UNIT=12 ,FORM= 'UNFORMATTED ' ,STATUS= 'OLD ' , 
1  ACCESS='DIRECT',RECL=LREC, 
1  IOSTAT=IERR,ERR=90) 

IF(IERR.NE.O)  GO  TO  90 
C 

READ(LUIN,REC=1)  NXP1,NYP1,NZP1 

NXY1=NXP1«NYP1 

IF(NXYI.GT. NXYIP)  GO  TO  92 

NREC=NXY1/LENREC 
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FILE:  FORMIZ   JCL      Al 


1SH0RT=NXY1-LENREC»NREC 

IF(ISHORT.GT.O)  NREC=NREC+1 
IF(ISHORT.EQ.O)  ISHORT=LENREC 

IF(DEBUG)  WRITEILUBUG,»)  NXP1,NYP1,NZP1,NXY1,NREC,ISH0RT 
C 

C   OPEN  ouipui  file:  FORMATTED,  SEQUENTIAL  and  of  FIXED  RECORD  LENGTH. 
C 

OPEN( UNIT=13 , FORM= ' FORMATTED ' ,STATUS=  *  OLD ' , 
1  ACCESS= 'SEQUENTIAL ' ,I0STAT=IERR,ERR=91 ) 
IF(IERR.NE.O)  GO  TO  91 
WRITE ( LUOUT, 100 )NXP1,NYP1,NZP1 
NR=1 

DO  £  K=1,NZP1 
DO  6  IC=1,3 
IMAX=LENREC 

DO  7  JR=1,NREC 

IF(JR.EQ.NREC)  IMAX=ISHORT 
NR=NR+1 

NRA=( JR-1)*LENREC 

READ(LUIN,REC=NR)  (AREC(NRA+I  ),I  =  1,IMAX) 
7   CO^mNUE 
C      Write  daia  "to  outpu't  file 

WRITE! LUOUT,101)(ARECl  II), 11=1, NXYl) 
6   COfrriNUE 
5   COhfTlNUE 

IFl DEBUG)  WRITE! LUBUG, 1900) 

1900  FORMAT (IX, 'NORMAL  STOP  IN  PROGRAM') 
GO  TO  999 

C  Error  "trapping.  .  . 

90  IF( DEBUG)  WRITE! LUBUG, 1901 ) 

1901  FORMAT! IX, 'ERROR  IN  OPENING  INPUT  FILE') 
GO  TO  999 

91  IF! DEBUG)  WRITE!  LUBUG, 1911) 

1911  FORMAT! IX, 'ERROR  IN  OPENING  OUTPUT  FILE') 
GO  TO  999 

92  IF! DEBUG)  WRITE! LUBUG, 1921 ) 

1921  FORMAT! IX, 'GRID  TOO  LARGE)  RESET  NXYIP  AND  RE-CREATE', 
1  /, IX, 'LOAD  MODULE'  ) 
999  CONTINUE 
C  CLOSE  files  ... 

CLOSE! LUIN,STATUS= 'KEEP'  ) 
CLOSE! LUOUT, STATUS= 'KEEP'  ) 
STOP 
END 
//60.FT12F001  OD  DSNsMSS.SlB-il .PH0ENICS.DF12,DISP  =  (0LD  ) 
//G0.FT13F001  DD  DSN=MSS.S1541.PH0ENICS.DF12C,DISP  =  !0LD ) 
/« 
// 
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